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THE ABSORPTION MICROWAVES EXCITED 


Abstract 


Direct absorption measurements were made upon excited atomic hydrogen 
discharge tube placed inside wave guide. Two ranges wave lengths the 
cm. region were investigated. The results can satisfactorily explained 
absorption due free electrons. Although the transition 


atomic hydrogen this region, selective absorption was observed. 


Introduction 


The energy levels the hydrogen atom for the second quantum state 
calculated the Dirac theory are shown Fig. Since transitions are 
allowed between levels the same principal number then according 
the selection rules for electric dipole radiation, transition between the 
and the levels permissible. This transition has wave number 


2 


9,914 mc. 
10,940 

+1034mc, 
shift 

Dirac. Lamb. 
(a) (b) 


Fic. Energy levels the hydrogen atom for the quantum state predicted Dirac 
and determined Lamb. 


0.365 which corresponds wave length 2.74 cm. Thus microwave 
energy wave length 2.74 cm. frequency 10,940 mc. should absorbed 
atoms the state, causing them jump the state. 
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Recent experiments performed Retherford and Lamb (11) indicate that 
the energy levels are shifted from the positions predicted the Dirac theory. 
This shift has been shown due interaction the electron with the 
radiation field (2) and resulting levels are shown Fig. This indicates 
that electromagnetic wave frequency 9914 mc. would absorbed. The 
results Retherford and Lamb show better agreement with the fine structure 
components the line found spectroscopically. Although most prob- 
able that the energy separation close 9914 mc., absorption measurements 
were also made the region 10,940 mc. 


The amount absorption that might reasonably expected was calculated. 
This required estimate the fraction the total number molecules 
the discharge which give rise atoms the state. Estimates this 
fraction for similar hydrogen discharge reported Carst and Ladenburg (4) 
and Haase (6) differ factor calculation based Carst and 
Ladenburg’s estimate, the more favorable the two, indicated that the ab- 
sorption might just the borderline possible measurement. The impor- 
tance direct measurement this energy separation justified attempt 
measure the absorption even though success was doubtful. 


The natural breadth the absorption line was calculated from the mean 
lives the energy levels from which arises. was found 100 mc. The 
Doppler broadening the line less than mc. which negligible compared 
with the natural breadth. Since the line breadth great, the conventional 
sweep frequency methods microwave spectroscopy cannot used. 
therefore necessary use fixed frequency methods. Any method producing 
excited hydrogen atoms large numbers will also give rise ions and free 
electrons. These ions and free electrons will absorb microwave energy resulting 
continuous background absorption. 


Early attempts other workers (3, detect microwave absorption 
excited hydrogen were unsuccessful. The development microwave tech- 
niques made possible more refined absorption experiments. Lamb and Reth- 
erford (8, obtained their results microwave method but did not meas- 
ure absorption directly. 


Apparatus 


The experiment required the production high concentration excited 
hydrogen atoms possible container through which microwave energy 
could passed. electrodeless discharge tube placed inside wave guide 
was used for this purpose. The amount energy absorbed was measured 
two ranges frequencies. Fig. block diagram the whole apparatus. 


Microwave energy was supplied klystron pulse modulated fre- 
quency c.p.s. with pulse about microseconds wide. The modulator 
supplied keying pulse for high speed oscilloscope which output levels 
were displayed. 


| 
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The main components the wave guide transmission line and measuring 
equipment block are shown Fig. The variable attenuator was calibrated 
the frequencies used and measured differences attenuation with 
accuracy considerably better than 0.1 db. small fraction the energy 


Photronic Cell 
Monitor for 
Discharge Tube 


Wave Guide 


Modulated Discharge Tube Termination 
Microwave Section Wave Guide 
Signal Equipment Wave Guide System 


Radio -frequenc 
Vacuum System 


Excit isch and 


Fic. diagram the apparatus. 


Oscilloscope Amplifier 


Variable Stub Fixed Flap Variable Standin 
Transformer (10 db) (calibrated) Detector 
Section 


Coaxial Cable 


Guide 
Coaxial-to- Variable Stub Tunable 
Transformer Wave Guide Impedance Crystal 

MONITOR BRANCH 
Signal 
Wavemeter 

Oscilloscope 


Fic. Block diagram wave guide transmission line and measuring equipment. 
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was fed through directional coupler monitor branch. The standing wave 
detector which followed the directional coupler was used conjunction with 
the variable attenuator for all standing wave measurements. 


Impedance Circular Impedance 


Transformer Wave Guide Transformer 


Mode 
Suppressor 


Section 
Discharge (Fig. Suppressor Wave Guide Section 


Fic. Discharge tube section wave guide. 


The components the next block Fig. containing the discharge tube 
are shown Fig. The size the circular wave guide surrounding the dis- 
charge tube was chosen that energy could transmitted the mode 
only. The rectangular guide was transformed the circular guide tapered 
section. 


External electrodes were used avoid much possible placing obstruc- 
tions the wave guide. These electrodes were insulated sections wave guide 
placed either end the discharge tube shown Fig. avoid leakage 
microwave energy the spaces insulating the electrodes, circular chokes 
and flanges were used. 


To Vacuum 


Flange Electrode System 


Discharge 
Tube 


Split Wave Circular 
Guide Wave 
Section Guide 


Fic. Discharge tube and external electrodes showing split wave guide section. 


connection the pumping system was brought out through small 
opening section split wave guide following the second electrode. Owing 
imperfections the circular wave guide, elliptical modes were produced. 
These were minimized much possible rotating the individual sections 
circular guide and the remaining effects were eliminated means mode 
suppressors placed either end the circular section containing the discharge 
tube. 


For measurements using the standing wave method, the wave guide system 
was terminated means solid plunger whose position could varied. 
For measurements the comparison method the wave guide was terminated 
the detector system shown Fig. 


From 
Standing A 
Wave 
Detector 
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The discharge tube shown Fig. was made mm. O.D. Pyrex 
tubing with tapered ends avoid reflection microwave energy. The con- 
nection the vacuum system was made through small diameter side tube 
mentioned above. The inside walls the discharge tube were coated with 
solution phosphoric acid before evacuation. This increased the concen- 
tration atomic hydrogen. 


Circular-to- 
Rectangular 
Wave Guide 

Taper 


Variable Stub 
Impedance 


Transformer 


Reflectionless 
Termination 


Directional 
Coupler 


Coaxial Cable 


From 
Mode 


Suppressor 
Coaxial to- Variable Stub Tunable 
Guide Impedance Crystal 

Transformer Transformer Detector 


Oscilloscope Amplifier 


Fic. Termination wave guide system for comparison method attenuation measurement. 


The discharge tube was excited radio-frequency oscillator capable 
delivering about 1.5 kw. mc. The voltage was stepped air-cored 
transformer before being applied the external electrodes. The radio-fre- 
quency current was measured with thermocouple meter. Operating current 
was from 250 ma. 350 ma. and the applied voltage was from 3000 4000 
v.r.m.s. 


The light output the discharge tube was monitored with red filter 
and photronic cell. The pressure the vacuum system ranging from 0.1 
0.4 mm. mercury was read thermocouple gauge. The hydrogen was 
produced electrolysis and dried over phosphorus pentoxide before admission 
the discharge tube. 


With the system described above, the discharge the tube extended uni- 
formly throughout the section between the electrodes for guide sections 
one meter length. 


Methods Measurement 


Two distinct methods were used measure the attenuation the plasma 
within the discharge tube, the standing wave method and the comparison 
method. 


i/ 
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Standing Wave Method 

This method based upon measurement the standing waves produced 
short circuit placed the end the wave guide beyond the section whose 
attenuation measured (5, 10). the plunger shifted the standing 
wave ratio will pass through maximum and minimum values. they are 


formula 


This method accurate only for fairly the input the 
attenuating section owing the experimental difficulties detecting small 
changes large standing wave ratios. The voltage standing wave ratio due 
the mismatch the input given 


10 log io db. 


denoted and ry;, respectively then the attenuation given the 


This formula used check that the match remains within certain limits. 


The attenuation was first measured with the discharge off. This was done 
for number frequencies each range. Then the measurements were re- 
peated with the discharge tube operating. The net attenuation due the plasma 
the discharge could then found. 


Comparison Method 


Attenuation was measured comparing the power entering the attenuating 
section with the power leaving the attenuating section. This method required 
that the whole system very well matched all points. 


Anomalous Dispersion Method 


attempt was made see the anomalous dispersion could detected 
the frequency where absorption was predicted. the vicinity absorp- 
tion frequency one might expect changes the dielectric constant. The value 
depends the dielectric constant and the positions the maxima and 
minima the standing wave pattern depend the value The effect 
anomalous dispersion might detected plotting the positions 
standing wave minima maxima against frequency. would show 
irregularity otherwise smooth curve. With the apparatus arranged for 
the standing wave method attenuation measurement, the plunger was left 
fixed position and the frequency changed. The positions the maxima were 
found for each series frequencies, with the discharge tube operating. 


Results 


Standing Wave Method 


The curves attenuation for the frequency ranges 9550 mc. 10,200 mc. 
and 10,750 me. 11,100 mc. are shown Figs. and Readings which were 
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Discharge 
mm. pressure 


Discharge 
pressure 


discharge 


9500 9600 9700 9800 9900 10,000 10,100 10,200 
Freq. 


Fic. measurements for the frequency range 9950 mc. 10,200 mc. with and 
without the discharge. 


Discharge 
-18 mm. pressure 


discharge 


Freq. mc. 


Fic. measurements for the frequency range 10,750 mc. 11,100 mc. with and 
without the discharge. 
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repeated agreed generally within 0.02 db. The net attenuation due the 
plasma found from these shown Figs. and 10. The deviations the 
points from the experimental curves were used calculate the probable error. 
This indicated the dotted lines. From these curves can concluded 
that for the frequency range 9550 mc. 10,200 mc. there selective ab- 
sorption greater than 0.06 db. which twice the probable error. For the 
frequency range 10,200 mc. 11,100 mc. there selective absorption 
greater than 0.14 db. 3.3%. 


mm. pressure 


~ 


Freq. mc. 


Fic. the plasma for the frequency range 9950 mc. 10,200 mc. the 
standing wave method. 


Freq. mc. 


Fic. 10. the plasma for the frequency range 10,750 11,100 mc. the 
standing wave method. 


Comparison Method 


The attenuation measured the comparison method for the lower fre- 
quency range shown Fig. which the dotted lines show the results 
the standing wave method the same range. The agreement between the two 
methods good and the accuracy about the same but the standing wave 
method very much easier apply since exact matching not necessary. 
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approximate calculation was made the amount background absorp- 
tion that would arise from the free electrons the discharge. This was based 
upon estimate the number free electrons per cubic centimeter 
obtained from value dielectric constant taken from the measurements 


Ain 
db. 


Freq, 
Fic. 11. the plasma for the frequency range 9950 mc. 10,200 mc. the 


comparison method. The limits for the attenuation the standing wave method are shown the 
dotted lines. 


used the search for anomalous dispersion. value for the number 
collisions made electron one second, was obtained from reported 
measurements hydrogen discharge tubes (1, 6). The absorption coefficient 
per centimeter (12) then given 


Position 
Maxima 
cm. 


10,095 10,050 9915 
Wavemeter Reading 


Fic. 12. Curve taken investigate anomalous dispersion. 


20 
18 
7 
16 
5 
4 
| 
| 3 
| iT] 
| 
| 
| 


CANADIAN JOURNAL PHYSICS. VOL. 


where the wave length the incident energy. This calculation showed that 
the background attenuation obtained may reasonably accounted for 
the absorption due free electrons. 


Anomalous Dispersion Method 


The curve showing position maxima the standing wave pattern plotted 
against frequency shown Fig. did not reveal any irregularities 
that could attributed anomalous dispersion. Small periodic departures 
from the smooth curve are attributed residual mismatch. The position 
these could changed moving the plunger. calculation the deviation 
from the smooth curve expected from quantum mechanical dispersion 
theory indicates that there only very slightly better chance observing 
the transition this method than attenuation measurement. 


Conclusion 


The estimate the success the experiment involved the population 
hydrogen atoms the state. The approximate value the population 
reported Carst and Ladenburg (4) now appears have been too 
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METEOROLOGICAL EFFECTS COSMIC RAY INTENSITY 
AND THE MESON 


Abstract 


Continuous recording the intensity the hard component cosmic rays 
has been carried out Ottawa (Geomagnetic Lat. 56.8) and Resolute (Geo- 
magnetic Lat. 83) the Canadian Arctic for over year. Two counter telescopes 
with different amounts lead shielding are used each station. The results 
have been analyzed for barometer effect four series observations each ex- 
tending over three months. When averaged, they show barometer coefficient 
per cm. mercury for mesons with momenta above 350 Mev./c and 
per cm. mercury for those with momenta above 630 Mev./c. 


calculation has been made the temperature and barometer 
was based the variations the probability decay mesons the atmos- 
phere, under varying atmospheric conditions. The calculations the temper- 
ature effect not yield results comparable with measurement; those the 
barometer effect show that considerable information can obtained about the 
meson spectrum from measurements the barometer coefficients. The results 
favor spectrum agreement with parts spectra measured other authors. 
They suggest the possibility second maximum the production spectrum 
mesons. 


Introduction 


The ‘‘barometer and the intensity cosmic 
rays sea level have played important role for many years (13) studies 
the nature cosmic rays. The temperature effect particular was shown 
Blackett (2) give proof the radioactive decay the meson. For some 
time both effects been considered consistent with current knowledge 
the production mesons the atmosphere and their decay and absorption 
they proceed sea level. more quantitative explanation, however, based 
this knowledge and actual atmospheric conditions seemed worth attempting. 


new calculation therefore has been made and measurements have been 
taken which were designed that their results could applied the calcu- 
lation. The calculation based the effect changing atmospheric conditions 
the probability decay mesons they pass through the atmosphere. 
actual mean atmosphere was used and alternative sea level spectral distri- 
butions were chosen compare with measurements. This approach differs 
considerably from that Janossy (8) who assumed isothermal atmosphere 
and chose simple mathematical expression represent the spectrum, from 
that earlier analyses (for instance Rossi (16)) wherein the effects absorp- 
tion the atmosphere and radioactive decay were treated separately. 


II. Experimental 


Description Apparatus 


Continuous recordings cosmic ray intensities have been taken two 
stations, Ottawa (Geomagnetic Lat. 56.8°N) and Resolute (Lat. 


received September 1950. 
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Long. 94°55’W. Geomagnetic Lat. 83°N) the Canadian Arctic. Both stations 
are near enough sea level (about 170 and ft. respectively) that consider- 
ation need given differences altitude. The Ottawa station has been 
recording almost continuously since October 1948 and the Resolute station 
since September 1949. Counter telescopes are used and, since September 1949, 
the equipment the two been identical. consists three counter 
trays arranged form two vertical telescopes two counter trays each, one 
tray being common both telescopes. The trays each contain two counters 
in. diameter and about in. active length. The trays are mounted one 
above another with in. lead horizontal sheets between the top and 
center tray (Nos. and and in. between Nos. and Double coin- 
cidences were recorded between trays and and between and There are 
therefore two telescopes, one with in. lead, the other with inches 
lead (plus negligible amount aluminum and brass) between trays. This 
arrangement was chosen give maximum amount information with 
minimum amount equipment, the space and power available Resolute 
were very limited. 

The bottom tray was placed floor level near the ground with in. lead 
under and in. each side. This shields effectively against soft com- 
ponent. The center tray and top tray were not shielded the sides, and the 
covering over the top tray and whole apparatus including the roof the 
building was equivalent gm. per wood, which was neglected 
comparison the lead shielding between counters. 


Previous September 1949, the arrangement counters Ottawa was 
somewhat different from that described above. For the first series (October 
1948 January 1949) two separate telescopes were used. Each had in. 
lead between counters and in. above the top tray. There was least in. 
lead the sides and the whole was placed inside old rifle range with con- 
crete walls and ceiling ft. thick. The second series (February June 1949) 
was similar the first except one telescope had in. lead between counters 
and the other had in. 

The results described here therefore come from four series measurements 
follows: 

Series 1—Ottawa, October 1948 January 1949 inclusive; shielding in. 

lead and ft. concrete; heavy shielding sides. 

Series 2A—Ottawa, February June 1949 inclusive; shielding 514 in. lead 

and ft. concrete; heavy shielding sides. 

Series but with in. lead plus concrete. 

Series 3A—Ottawa, September December 1949 inclusive; in. lead 

between counters and light roof. 

Series 3B—Same but with in. lead between counters; bottom 

counter shielded sides. 

Series same period and arrangement 3A. 

Series 4B—Resolute, same period and arrangement 3B. 
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The minimum momentum accepted the shielding was estimated from the 
curves published Montgomery (14). 


Telescopes with only two counter trays and side shielding not give high 
discrimination against side showers. This important because the known 
high barometer effect for showers. However, calculations and tests showed that 
the experimental arrangement used here the number coincidences due 
side showers was less than and the error due these, even though 
systematic, would not affect the conclusions. Therefore, corrections were 


made for side showers. 


The method recording was developed especially for long periods con- 
tinuous operation and very accurate. Multiple binary scalers were used and 
each time scaler repeated its cycle was marked multiple pen recorder 
which time marks were included. The details the circuits and methods 
recording will published separately. 


Methods Analysis 


the same records were used study anomalies intensity due such 
phenomena solar flares (15) the records were averaged over short intervals, 
usually four hours. save computation and make best use statistical 
accuracy counting rates some selection was made the data used. 


For analyses with barometric records, periods eight hours more were 
chosen during which the barometric pressure was reasonably steady, that is, 
did not change more than about millibars during the period. This repre- 
sented less than half the total time. The mean cosmic ray intensity during these 
periods was compared with the pressure standard statistical methods. The 
selection data from periods steady pressure seems justified because points 
taken during rise fall the barometer will give values near the mean where 
there are large number points any case. The selected data are dis- 
tributed more uniformly over the range pressure, and the amount com- 
putation greatly reduced. Occasionally shorter periods (less than eight hours) 
were included order take advantage very low pressures which not 
rule last for long periods, and one series measurements where hr. 
averages had been worked out these were used directly, again selecting those 
where the barometer was reasonably. steady during the period. 


the correlation with barometric pressure temperature corrections were 
made. Forbush (7) has shown that correction for ground temperature does 
not increase the reliability the data. The objective this measurement was 
get barometer coefficient that could compared with the theoretical 
calculations Section III. this the measured barometer coefficients 
should independent temperature. Though temperature changes through- 
out the atmosphere obviously affect the cosmic ray intensity sea level, 
straight correlation with pressure should give the true barometer coefficient 
the relation between temperature and pressure random nature. 
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There well known inverse correlation between pressure any specific 
height above sea level and the mean temperature between that height and sea 
level. However, the analyses the next section the behavior mesons 
studied actual atmosphere specified terms pressure and density 
rather than terms height. rough test was made look for correlation 
between ground pressure and temperature any pressure level the atmos- 
phere. Data were obtained from the radiosonde meteorological station Sault 
Ste. Marie, which the nearest radiosonde station Ottawa with 
weather. Using the periods steady pressure selected for the cosmic ray mea- 
surements Series the excursions from the mean seasonal temperature 
various pressure levels 200 mb. were plotted against ground pressure. 
The resulting plot showed inspection such random scattering points 
that the correlation seemed negligible. The barometer coefficient obtained 
straight correlation with ground pressure therefore probably the best 
compare with the calculation Section III even though some seasonal tem- 
perature effects may not have been eliminated. 


The analysis used here was planned for comparison with the calculations 
Section III and not particularly for the correction cosmic ray data. Therein 
differs from that Duperier (6) who finds very good correlation between 
the intensity the hard component sea level and combination three 
variables, the barometric pressure ground level, the height the effective 
meson production layer (the 100 mb. layer), and the mean temperature 
levels between 100 and 200 mb. the atmosphere. Duperier’s analysis would 
give more precise way correct cosmic ray data where information the 
upper atmosphere available. 


Results 


The measured barometer coefficients are shown Fig. and listed Table 
with corresponding correlation coefficients. The shielding was such that the 
points are distributed about two groups, one average momentum 351 Mev./c, 


TABLE 
Barometer Momentum 
Series coefficient limit, 
%/cm. Mev./c. 
305 
374 
420 
Mean 351 
-1.88 604 
682 


Mean 630 
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and the other 630 Mev./c. The means are shown these two points Fig. 
The distribution the points about the means gives indication the 
accuracy. This not high, but the results suggest clearly that the magnitude 
the barometer coefficient increasing with momentum over the range mea- 
sured. This feature consistent each series observations where two thick- 
nesses absorber were used simultaneously, and will seen Section III 
have important bearing the shape the meson spectrum low 
energies. 

This upward slope toward low momenta not entirely agreement with 
other observers where different thicknesses absorber were used. Most mea- 
surements recorded with two thicknesses absorber were taken with ioni- 
zation chambers. Though measurements the barometer effect taken with 
ionization chambers and with counters should vary the same way with mo- 
mentum, they are not quantitatively comparable. The counter telescope point- 
ing vertically includes only limited cone particles and counts each one with 
equal emphasis. chambers have directional selection and over- 
emphasize the particles that just penetrate the width the chamber. Some 
values obtained with ionization chambers are: 


Compton and Turner (5) cm. lead 
Messerschmidt (10) cm. lead 


BAROMETER COEFFICIENT, cm. 


500 1000 1500 2000 


MOMENTUM, 


Fic. The integral barometer coefficients. The points represent measured values and the 
curves computed values. Points each series measurements are marked alike. represents 
the means groups points. Curve corresponds spectrum for whichy Curve the 
Wilson spectrum, Curve spectrum for which and Curve was selected fit the 


measured points. 
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Values with counter telescopes are: 
Duperier (6) cm. lead 
Barnothy and Forro (1) cm. lead 
Kolhorster (9) 140 wood 
Smith and Gast (12) 20.5 cm. lead 
(reported 30. cm. lead 
Millican and Loughridge) 


Duperier has pointed out the great variation results obtained when 
simple correlation made between intensity and pressure; this was also noted 
here the analysis the results Series was made monthly periods 
well for the whole period four months. The monthly values the baro- 
meter coefficients are —2.52, and and the correlation 
coefficient —.66, and respectively. simple correlation this 
type value obviously must carried out over long period. 


not certain whether any significance should attached differences 
between measurements Ottawa and Resolute. While the mean climate 
very different from the point view comfort, will not have much effect 
the probability meson decay the atmosphere. The sea level intensity 
essentially the same Ottawa and Resolute, both are above the 
the latitude effect. Significant differences barometer coefficients lower 
latitudes might expected. 


Theoretical Calculations 


General Considerations 


Barometer and temperature coefficients may calculated from knowledge 
the sea level spectrum and the behavior mesons their passage through 
the atmosphere. Rossi (17) has given reasons for assuming that the flux 
ionizing particles the hard component (under cm. lead shielding) sea 
level consists almost entirely mesons. The small fraction heavier particles 


and few electrons that may come from showers the lead can 


This has been accepted the following calculations even though small frac- 
tion penetrating particles with large meteorological effects (such are 
found for penetrating showers) would affect the results. The process the 
formation mesons (the decay meson meson) also not con- 
sidered relevant view the very short life time the meson. 


assumed that mesons reaching sea level with momentum have come 
the average from pressure level the atmosphere, the level being 
measured its pressure grams per square centimeter. The probability 
their decay between production level and sea level can calculated from the 
well known radioactive decay equation, the form 
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where represents the number particles any point with momentum 
the mass and the mean life rest, and element length path 
the meson. For mesons moving vertically downwards may written dp/p, 
where the density the atmosphere. one assumes that the momentum 
loss per unit mass air traversed constant over any range momenta from 
production level sea level (momenta below minimum ionization are excluded 
with the soft component) may written p), where the 
momentum sea level and the momentum loss per gram per square centi- 
meter. The ratio intensity sea level that production level 


Num 
The the denominator the exponent will normalize all units that mass 
measured units Mev./c?, momentum Mev./c, pressure grams per 
square centimeter, density grams per cubic centimeter, and seconds. 
The path the particle assumed rectilinear, nearly rectilinear 
that scattering has negligible effect its range. The occasional collision 
where the energy loss very large also neglected. 


change atmospheric conditions will result change the original 
momentum meson which reaches sea level with momentum Pg. pro- 
duction spectrum the inverse power law form assumed and the intensity 
sea level can then represented the equation 


where normalizing constant and the energy production. This 
equation cannot integrated directly unless some simple relation between 
pressure and density exists. also shows the difficulty calculating separate 
barometer and temperature coefficient. This was done however selecting 
standard atmosphere, then second atmosphere differing from the first 
temperature only, and third differing from the first pressure only. 


The standard atmosphere chosen was the mean (from radiosonde measure- 
ments Sault Ste. Marie, Michigan) the atmospheric conditions during 
periods steady barometer selected Series for the analysis cosmic ray 
results. The temperature-altered atmosphere was calculated from the standard 
atmosphere varying the temperature amount was chosen 
for AT) each pressure level the atmosphere. This not real atmosphere, 
temperature changes never extend uniformly all the way through the 


atmosphere. 


The pressure-altered atmosphere was calculated changing the pressure 
amount sea level and proportional amount other levels. 
The density was adjusted that the temperature any height above sea level 
remained constant. This also may not represent pressure change found 
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practice but temperature changes throughout the atmosphere vary 
random way with the total weight atmosphere this pressure-altered atmos- 
phere compared with the standard should give barometer coefficient com- 
parable with the measurements described Section 


The Temperature Coefficient 

Equation (1) was integrated numerically for the three atmospheres, each 
with series values for Pg. The temperature coefficient defined 
and usually quoted per cent per degree Celsius. the value the expon- 
ential part the equation after integration represented for the standard 
atmosphere and for the temperature-altered atmosphere, the temperature 


coefficient may written 


IAT SAT 


and refer the energies production mesons which reach sea level 
with momentum the total mass air through which the mesons pass 
will the same for both the standard atmosphere and temperature-altered 
atmosphere Therefore, the temperature coefficient 
independent the production spectrum. Actual measuring apparatus usually 
accepts particles with all momenta above the minimum fixed the shielding. 
Therefore, coefficient that required. This was obtained 
further numerical integration based the calculated differential coefficients 
and the sea level spectrum. For the temperature effect the sea level spectrum 
Wilson presented Rossi (17) was chosen. Fig. shows the differential 
temperature coefficient calculated from Equation (2) and the integral curve. 
The coefficients are negative, would expected..As mentioned above, the 
temperature-altered atmosphere does not represent actual conditions. Apart 
from repeating all the numerical integrations for every change atmospheric 
conditions there way compare calculations with measurements tem- 
perature effects. Duperier (6) has shown partial correlation between sea level 


INTEGRAL 
DIFFERENTIAL 
- 
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Fic. The differential and integral temperature coefficients based uniform change 
temperature throughout the atmosphere. 


: 
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intensity and temperature levels between 100 and 200 mb. the atmos- 
phere. calculation could have been made for temperature change that 
level only, but would give small negative coefficient rather than the positive 
value found Duperier. 


The Barometer Coefficient 


the exponential part Equation (1) again represented for the 
standard atmosphere and for the pressure-altered atmosphere, the baro- 
meter coefficient may represented the equation 


this case not equal Er. Therefore, the spectral term important. 


cm.Hg 


BAROMETER COEFFICIENT, 


1000 2000 3000 4000 
MOMENTUM, 


Fic. The calculated differential barometer coefficients for different values Curve 
represents spectrum selected fit the measured points Fig. 


Fig. shows the calculated differential barometer coefficients for 
and for spectrum (marked which gives good fit with the measured 
coefficients. mentioned the case the temperature coefficient integral 
coefficient required for comparison with measurements. The integral baro- 
meter coefficients were calculated further numerical integrations over 
selected sea level spectra using the differential coefficients calculated from 
Equation (3). The results are shown the curves Fig. 1.* For Curve the 


The units used for pressure the calculation barometer coefficients were grams per square 
centimeter. The coefficients are usually quoted per centimeter mercury. Therefore, this unit 
has been used displaying results. Similarly pressure levels the atmosphere are usually 
measured millibars. Hence that unit used referring production levels mesons. 


J 
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chosen sea level spectrum that obtained when The calculations for 
Curve were based the sea level spectrum measured Wilson (as pub- 
lished Rossi (17)). Curve calculated from spectrum chosen give 
good fit with the measured points and corresponds Curve Fig. For 
Curve course applies the production energy spectrum while the 
final integration required calculate the curves Fig. based the sea 
level momentum spectrum. The data from the integration Equation (1) can 
used derive one spectrum from the other. 


IV. Discussion 

General Remarks 

Interesting features the results these calculations are the values the 
differential barometer coefficient for and the way which this coefficient 
varies with the production spectrum showed the same intensity all 
energies the value the differential barometer coefficient would positive 
low momenta instead negative. This can explained the fact that, for 
low values the change mean momentum which results from change 
barometric pressure overshadows the change the height the effective 
layer formation. The sensitivity the integral barometer coefficient the 
slope the production spectrum represented the value shown 
the difference between Curves and Fig. 


Accuracy Calculations 

The values used for the mass and mean life rest for the meson were 215 
times the electron mass and 2.10 sec. respectively. the calculated 
barometer coefficients represent differences two integrations Equation 
(1), possible errors these would make little difference the results. 


When the calculations were started was assumed that would suffi- 
ciently accurate treat constant having value slightly above that for 
minimum ionization. rough mean between 500 and 1000 Mev./c was chosen, 
making 1.95 Mev./c per gm. per Early results showed however that 
the part the meson spectrum between 2000 and 10000 Mev./c. played 
dominant part fixing the general level the integral barometer coefficient 
curve. The above value was too low for that region. Therefore, correction 
was made recalculating some points using values line with published 
data for the mean momentum for the value concerned. 


The value the mean level production, was also assumed first 
constant and fixed 100 mb. calculating the production spectrum from 
the sea level spectrum, Sands (18) has pointed out that the mean level pro- 
duction must vary with momentum, particularly the low end the spec- 
trum. assumed production law the atmosphere the form 
grams per square centimeter and calculated the effective level production 
for mesons different range. This method was not used here, but the differ- 
ential barometer coefficients were calculated for assumed production levels 
200 and 300 mb., well 100 mb., establish the effect varying heights 


| 
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the results. Finally production height function was adopted wherein 
mesons with sea level momentum 200 Mev./c (the lowest considered 
these calculations) were assumed come from mean level 150 mb. The 
height production was assumed slope upwards 100 mb. for mesons with 
sea level momentum 1000 Mev./c and remain constant for those 
higher momentum. This height function represents slightly higher levels than 
that used Sands but the calculations show that, though the height function 
has considerable influence the exact shape the derived production spec- 
trum, the differences would not affect the general conclusions that can 
drawn from the measurements described here. 


The derivation Equation (1) such that applies only mesons moving 
vertically downwards. The counter telescopes used included quite wide 
solid angle (in the case the shortest much 10° 50° half angles). The 
effect this can shown small rederivation Equation (1) for 
mesons whose path makes angle with the vertical. first approximation 


cos cos 


and therefore does not change. The ratio changes very little. The largest 
error occurs replacing This infers linear relation between 
the range lead, and momentum right down low momenta. also applies 
only arrangements wherein the shielding horizontal sheets, which was 
the case the measurements described here. With cylindrical spherical 
symmetry shielding, the directional effects would more pronounced. 
They cannot however very great because the case Series analysis 
was made the recordings counting rates from the bottom counter telescope 
which had in. lead top pile in. wide but only enough the 
sides remove the soft component. The mean value the barometer co- 
efficient was per centimeter mercury. rough calculation one 
value shows that, with the counter telescopes used, the error due 
solid angle would direction that would tend modify the conclusion 
(see Section about second maximum the production spectrum mesons 
but not great enough allow its withdrawal. exact computation in- 
cluding directional effects would involve further integrations and assumptions 
about differential spectrum different angles. 


The numerical integrations were carried out using trapezoids maximum 
width give reasonable accuracy withoyt excessive computation, and certain 
amount interpolation was used. Errors due the methods used however 
could hardly represent more than about 0.1% per centimeter mercury the 
position the curves Fig. and would not affect the conclusions. 


The Meson Spectrum 

Comparing the points Fig. with Curve for which seen 
that this curve passes about midway between the means the two groups 
measured points. The points, however, indicate downward negative slope 
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with increasing momentum, while Curve slopes upward. The accuracy such 
that the significance this downward slope questionable but the conclusions 
that result from accepting are interesting and have some support other 
experiments. 


simple power law spectrum will not give integral barometer coefficient 
curve with this negative slope low momenta. Rossi 1948 reviewed measure- 
ments the sea level meson spectrum and showed both differential mo- 
mentum and range spectrum (Figs. and ref. (17)). gave reasons which 
indicate that the measured points low momenta underestimate the intensity 
and showed solid curve representing the most likely sea level spectrum. 
Translating this production spectrum using the present calculations, 
spectrum obtained which gives almost constant value little less 
than low momenta and rising about the higher momenta values. 
This agreement with similar calculations made Sands (18) who gives 
average value 2.5 for the corresponding part his spectrum.* 


the measured points suggested slope that not agreement with 
simple power law production other spectra were tried. Curve Fig. 
based smooth curve drawn through Wilson’s measured points shown 
Rossi (Fig. ref. (17)). This gives the desired slope low momenta but 
general level well above the measured points. Curve was based 
spectrum which, momenta below 3000 Mev./c (at production levels), iden- 
tical that derived from Wilson’s measurements. However above 3000 
Mev./c has slope corresponding 3.5. This curve gives reasonably 
good fit with the measured points, and the increased slope high momenta 
agreement with recent measurement the sea level spectrum Caro, 
Parry, and Rathgeber (4). When transposed production levels their curve 
would indicate value close 3.5 above 3000 Mev./c, but lower mo- 
menta more accord with Rossi’s sea level spectrum. 


The production spectra corresponding Curves and are shown 
Fig. Spectrum suggests maximum the intensity about 2000 Mev./c. 
The intensity much lower momenta must much higher than the value 
this point account for the increase total intensity with height the 
atmosphere and consistent with measurements high altitudes (3, 11). 
Therefore, the suggested maximum 2000 Mev./c must second maximum 
the production spectrum. The present results are not considered sufficiently 
significant offer more than suggestion such second maximum but there 
some additional evidence its favor. The measurements Miller, Hender- 
son, Potter, Todd, and Wotring (11) altitude 3.4 km. show the ex- 
pected maximum low momenta, and their curve for positive particles 
under cm. lead there the suggestion asecond maximum slightly below 
1000 Mev./c. These authors not consider this maximum significant but 


Sea level measurements not give any information about the very low momentum end the 
meson production spectrum because, unless the mesons have enough energy reach sea level from 
height where appreciable number are produced, the contribution the sea level intensity 
negligible. the present calculations the lowest sea level momentum considered (200 Mev./c) 
represents momentum about 1900 Mev./c mean production level 150 mb. 


_ 
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reduction spectrum the expected spectrum 3.4 km. shows maximum 
very roughly agreement with curve through the observed points Miller 
al. The broken curve Fig. has therefore been drawn qualitative 
representation production spectrum that might expected momenta 
values below those included spectrum 


INTENSITY 


RELATIVE 


MOMENTUM, 


Fic. The production momentum spectra mesons. Curve corresponds constant 
value Curve the Wilson spectrum, and Curve the spectrum suggested measure- 
ments the barometer coefficient. 


Summary Conclusions 


The simple theory derived here gives new physical explanation the 
meteorological changes the sea level intensity the hard component 
cosmic rays. This based current knowledge the production mesons 
the atmosphere and their loss energy and probability decay passing 
through the atmosphere. 


Calculations the effect are necessarily based hypo- 
thetical changes the atmosphere and, therefore, not give results that can 
compared with measurements. They indicate, however, general agreement 
with the expected behavior. 
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Calculations the effect have been made which can com- 
pared with measurements. These show that the slope the meson spectrum 
production levels the atmosphere has dominant influence the baro- 
meter coefficient. The measured results therefore indicate some features the 
production spectrum that are agreement with parts spectra measured 
other observers. momenta below 1000 Mev./c there good agreement with 
the measurements the sea level spectra Wilson, but higher momenta 
the results show better agreement with the measurements Caro, Parry, and 
Rathgeber. 


Though the accuracy not great the production spectrum which best fits 
the barometer coefficient measurements suggests that maximum may exist 
about 2000 Mev./c, which must second maximum the total produc- 
tion spectrum. (The first much lower momenta.) This second maximum 
also suggested recent measurement the spectrum 3.4 km. Miller, 
Henderson, Potter, Todd, and Wotring. 
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PHOTOELECTRIC METEOR 


Abstract 


The photoelectric meteor detector consists type 931-A photomultiplier 
tubes arranged cover the entire sky. Signals from the photocells are re- 
peatedly sampled time sequence and displayed cathode ray tube. Photo- 
electric meteor signals, radar meteor echoes, and marker light operated 
visual observers are photographed continuous-film camera. 


Preliminary tests indicated that the average sensitivity the photoelectric 
detector meteors was about one visual magnitude inferior that the human 
eye. Approximate meteor positions and qualitative light intensities could 
obtained from the photoelectric record. Timing comparisons between the 
photoelectric and the radar records could made 0.05 sec. Methods im- 
proving the sensitivity and other operational characteristics are discussed. 


Introduction 


Since the summer 1947 the National Research Council and the Dominion 
Observatory have been engaged joint program meteor observing, em- 
ploying radar and radio equipment, visual observers, and both direct and 
spectrographic cameras for meteor photography (1, 2,3). More complete details 
the projects are available these earlier papers. During the observing 
periods the visual observers scan the sky and plot maps the visual meteor 
paths, while the same time record meteor echoes obtained 
from the radar system. When visual observer sees meteor presses button 
which flashes marker light the radar display. This marker the radar 
film has potential timing accuracy 0.1 sec. but, owing the reaction times 
the observers, the mean relative timing accuracy the order 1/2 
sec., and the individual error can considerably greater. Satisfactory cor- 
relation for most purposes achieved this technique the extent that 
definite coincidences between the brighter meteors and the radar echoes can 
well established, and can further confirmed comparing the visually 
plotted meteor positions with the positions obtained from three-station tri- 
angulation the radar echo ranges (4). 


However, there evidence that the radar echo path length shorter than 
the visual path general, and that the radar path segment may tend occur 
towards the latter portion the plotted visual path. Before definite con- 
clusions can drawn would desirable time the onset the meteor 
luminosity relative the radar echo with accuracy least 0.1 sec. and 
preferably 0.01 sec. The use photocells detect the meteor light has been 

suggested. the photocells were prove satisfactory this respect was 
expected that they might valuable other phases the work, for example, 
determine the general luminosity curve the meteor, particular time 
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the appearance bursts, flares, and give rough plot the meteor 
position and even the meteor path the sky. 


The photocell system described here has been moderately successful 
meeting the general specifications. During the trial periods observation 
several modifications have been suggested for improving some features the 
design, and these will incorporated future operations. However, the 
project has been temporarily suspended, was thought desirable present 
interim report outlining the progress date. 


Apparatus 


Each photocell unit consists 931-A photomultiplier tube 
associated battery power supply. Rectangular apertures have been cut each 
photocell box admit sky light over approximately 30° azimuth and 36° 
elevation. cover the entire sky these photocell units were arranged 
circular table shown Fig. The amplifiers and gating units can 
seen the cabinet supporting the photocell table. 


The block diagram, Fig. illustrates the general principles the system. 
essence, the information obtained the photocell units sampled time 
sequence the electronic circuits such manner that can readily 
correlated with the radar record. The initiating trigger pulse from the radar 
transmitter passes directly the range sweep circuit. Identical linear sweeps, 
260 km. long, are applied horizontally two cathode ray tubes located side 
side. The tubes are photographed camera which mm. film moves 
vertical plane continuously past the open lens in. per min. Seconds 
pulses appear each sweep facilitate the precise timing the phenomena 
appearing the two C.R.T.’s. Twenty-kilometer range markers are impressed 
the sweeps well, which serve the conventional purpose the radar dis- 
play, but the photocell display they are used merely gate markers 
identify the outputs the photocell units. 


The trigger pulse from the radar transmitter applied amplifier unit 
No. through electronic delay km. (267 sec.) because the range 
interval from km. reserved for the first row seconds pulses. This 
delayed trigger then opens No. electronic gate circuit for interval 
km. (67 sec.), allowing the light signal from No. photocell unit pass 
the photocell display C.R.T. soon the gate No. unit has closed, 
trigger pulse sent No. unit open the gate this unit for km., and 
the 19th unit. practice, additional delay inserted just prior 
the 19th unit avoid interfering with the second row seconds pulses 
appearing between 220 and 240 km. the range sweep. The gated outputs 
all the photocell and amplifier units are then combined and applied 
intensity-modulation the photocell C.R.T. display. Thus, the output from 
particular photocell will appear only the interval the range trace 
selected for the gating system. The signal from given photocell 
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sampled for sec. every 1/120th second (the transmitter pulse repe- 
tition frequency approximately 120 c.p.s.) but successive sweeps overlap 
the film that the light intensity record appears continuous. The practical 
limit time resolution about 0.05 sec. with this recording technique. 


Fig. shows some selected records the radar, visual, and photocell obser- 
vations. The radar record the upper one each photograph, the black 
horizontal lines being the km. range markers and the white dots the second 


photoelectric meteor detector. Nineteen photocells, covering the entire sky, are 
arranged the circular table. The cabinet contains the amplifier and gate units. 
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pulses. The diamond-shaped markers between the radar and the photocell 
records are impressed the film the visual observers. The photocell gate 
numbers are indicated the side the photocell record, shown the lower 
portion the photographs. radar echo, photocell signal, clearly dis- 
tinguishable the records virtue its size, shape, intensity, and duration, 
comparison with the spots and streaks random noise impulses present 
both the radar and photocell systems. Photocell signals due rotating 
airport beacon, located two miles west-south-west the observing post, appear 
sequence gates No. 12, 13, and intervals 9.5 sec. 


PHOTOCELL AMPLIFIERS GATES TRIGGER PULSE 


Fic. Block diagram the photoelectric meteor detector. 


The transient response photocell and amplifier was the order 
millisecond, which was considerably better than required. Condenser-coupling 
was employed the amplifiers, using long time constants. d-c. amplifier 
had been employed the background level, due mainly starlight, would have 
been high comparison with meteor light signal that means bucking 
out this component would have been needed. Furthermore, 
continuous corrections would have been necessary allow for slow changes 
the ambient light background. For short-duration light impulses approxi- 
mate reproduction the light intensity was obtained, subject nonlinearity 
and limiting effects the amplifiers and the C.R.T., but for longer duration 
light signals the order second more only the rapid changes intensity 
would show, owing the differentiating characteristics the amplifier. 


The battery photocells and the cabinet associated amplifiers and gate 
circuits are placed the open, well away from buildings obstructions. 
canvas fence protects the apparatus from very low angle illumination. Apart 
from the provision 110-v. 60-cycle power for the electronic equipment the 
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only connections between the outdoor photocell assembly and the display 
located indoors are two coaxial cables, one supply the trigger pulse the 
gating circuits and the other transmit the combined signal outputs the 


remote display. 


Observations 


The photocell system was tested intermittently during the autumn 1949, 
and the first operational run was made during the observations the Geminid 
shower, December 12-14, 1949. Three radar stations and one continuous-wave 
Doppler station were operation. During the hours darkness, direct meteor 
cameras and spectrographs were used and crew three four experienced 
observers plotted visual meteors under Dr. Millman’s direction. Examination 
the radar-photocell film records revealed appreciable number apparent 
time coincidences between radar meteor echoes and the photocell signals. How- 
ever, the present discussion will confined chiefly those meteors for which 
coincident visual observation was available, because the plotted position and 
path the meteor and its estimated visual magnitude are necessary anal- 
yzing the value the photocell observing technique. Twenty-two triple coin- 
cidences between radar, visual, and photocell observations were selected and 
are summarized Table Four meteors for which radar echo was recorded 


have also been included. 


The sky coverage the photocell system shown Fig. The zenith 
the center and the outer boundary approximately the horizon. The 
correlation between the visually plotted paths and the photocell records 
reasonably good allowance made for the observational errors. The average 
error the visual plotting the both azimuth and elevation 
but individual plot may have considerably greater error. The boundaries 
the individual photocell sectors are not sharp but have wide penumbra, 
and may furthermore error few degrees. 


The grouping the majority the meteors the western half the sky 
not accidental, nor does imply that the eastern photocells were operating 
inefficiently. general rule the average meteor will detected the radar 
system only when its elongation from the radiant, near 90°, though 
bright meteor may often detected any elongation. may shown (1) 
that this condition tends limit the radar observations those meteors 
occurring zone across the sky opposite the radiant. Between the hours 
1930 E.S.T. and 0100 E.S.T. the true bearing the Geminid radiant lay 
the eastern sector. Thus, most the triple coincidences between radar, photo- 
cell, and visual meteors should occur the western part the sky, since most 


the observed meteors were Geminids. 


The first column Table gives the mean visual magnitude the meteor 
seen the photocell system, and the second column lists the mean visual 
magnitudes all visual meteors recorded during the selected observing periods. 


d 


116 CANADIAN JOURNAL OF PHYSICS. VOL. 29 
TABLE 
= 1 =A { | } 
| | | (a) } plot 
| | | | 
Dec. 12, 1949 
Dec. 13, 1949 
Dec. 14-15, 1949 
TABLE 
Mean magnitude Mean magnitude 
Observing period photocell all visual 
meteors meteors | 8 
Dec. 12, 1949 2.3 1.2 
1937—2255 E.S.T. 
Dec. 13, 1949 1.5 3.8 
1928—2102 E.S.T. 
Dec. 13-14, 1949 1.9 
All periods 0.0 2.1 2.1 
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The sensitivity the photocell system during the December 12th’ period 
appears have been about magnitude less than that the eye. The data 
are weighted favor the visual observations because all these have been 
averaged, including the faintest meteors, but the photocell data include only 
the brighter ones for which identification was positive. The performance 
the photocells has dropped considerably the next two nights, partly owing 


Fic. coverage chart. This stereographic projection; the radii the inner, 
middle, and outer circles are respectively 18°, 54°, and coverage areas are indicated 
circled numbers. Solid arrows—plots photocell-visual-radar meteors. Dotted arrows—plots 
photocell-visual meteors. 


poorer observing conditions indicated similar, though smaller depre- 
ciation the visual magnitudes. The photocell battery voltages deteriorated 
with time and the cold weather and were below normal the last two nights, 
which would account for the poorer relative performance the photocell 


system. 
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While the sensitivity limit the photocells appeared about the second 
magnitude large number visual meteors the range zero second magni- 
tude were not definitely recorded the photocells. the other hand the 
photocells did record smaller number probable meteors which were not 
visually observed, although coincident radar echoes were sometimes noticed. 
similar effect has been remarked the correlation between radar echoes and 
visual meteors when the fainter magnitudes are considered. The radar-visual 
relationship not yet fully understood; often bright visual meteor yields 
weak meteor echo and vice versa. The principles production and detection 
radio echo from the ionized gas column left the meteor differ consid- 
erably from the production and detection visible light, but the main differ- 
ences between the photocell and the eye lie relative sensitivity and spectral 
response. The spectral response curve the 931-A photocell has maximum 
4000 and down 50% 3200 and 5400 which suggests that 
these photocells may discriminate against the slow meteors red orange 
color. The sensitivity the photocell used this application limited 
the variations starlight, light the night sky, aurora, and other ambient 
light, and not inherent noise the photocells amplifiers. The noise level 
output the photocell drops proportion size the aper- 
ture, and fainter meteor magnitudes could detected against the star 
background simply reducing the aperture. Adding more photocells with 
smaller apertures the system would not only improve the sensitivity but 
would also increase the accuracy position plotting, but there economic 
limit this procedure (the amplifier and gate units alone contained 140 vacuum 
tubes, exclusive the number the C.R.T. display unit). Alternatively, one 
could employ the same number photocells cover smaller portion 
the sky. 

Owing instrumental defect the gate circuits, which has since been 
rectified, the occurrence strong light signal one the photocell channels 
caused the sensitivity several the adjacent channels depressed 
momentarily, or, some instances, introduce fainter spurious signals 
other channels. More the meteors listed Table might have been recorded 
two more photocells were not for this effect. 


few detailed comments follow some the meteors listed Table 
and the photographs Fig. Reference should made Table and 
Fig. 


Meteors No. and 


The visual plots for these meteors lie adjacent the chart areas where 
they were seen the photocells. The visual plots may have been error 
few degrees, the coverage Photocell No. may have differed 
appreciably from the nominal outline shown Fig. The time dis- 
crepancies the case Meteor No. leave some room for doubt, though 
not unusual for the radar echo occur seconds after the passage 
the meteor, and the visual observers were occasionally late marking. 
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Meteor No. (Fig. 3a) 
The light signal from Photocell No. was coincident with the radar echo 
and endured for 1.2 sec., with indication flares toward the end. 


Meteor No. 
This was low the horizon that the radar echo, there were one, would 
have occurred greater range than the maximum 260 km. appearing 
the record. The photocell signal was strong enough suggest that the 
visual observer may have underestimated the meteor brilliance. 


The signal from Photocell No. has lasted 0.3 sec. The noise level Photo- 
cell No. has been set too high. 


Meteor No. 15. (Fig. 3c) 
The radar echo occured about sec. later than photocell signal. 


Meteor No. 17. (Fig. 3d) 
The photocells did not register till the meteor approached 90°, appre- 
ciably later than the visual observation the initial appearance the 
radar echo. 


Meteor No. 18. 
remark the recorder’s sheet indicates that the recorded visual time 


was several seconds late. 


Meteor No. 21. (Fig. 3e) 
The photocells have detected the luminosity later part the meteor 
path than have the visual observers. This confirmed the position 
the meteor plot Fig. Signals other gates are spurious. 


Meteor No. 23. (Fig. 3f) 
The radar echo was delayed 6.5 sec. after the visual and photocell times. 
This characteristic the echoes from bright meteors near the radiant 
27°-32° for this meteor), which are often several seconds late 
appearance, and may show multiple echoes characteristic (1)), and 
long duration; this case the echo lasted for 110 sec. 


Meteor No. 
Another example delay the appearance the radar echo when the 
observed meteor path close the radiant. 


will remarked from Table and the Chart, Fig. that Meteors Nos. 
,17, 20, 21, 23, 24, and all appear have registered the 
photocells toward the end the visual path, even beyond the plotted paths. 
The data for Meteor Nos. and suggest that the photocell detected these 
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meteors earlier than the visual observers. the majority cases the evidence 
indicates that the meteor luminosity increases end the visual 
path. 

Conclusions 


this stage the development the photoelectric meteor detector 
attempt will made interpret the observations exhaustively because 
imperfections the instrumentation. Several minor defects were found 
the electrical circuits, some which have since been remedied. The sky 
coverage each photocell was known only approximately. The addition 
long cone shield front each aperture would not only aid locating the 
field view more accurately but would greatly reduce the penumbra the 
edges the field, caused the present position the aperture couple 
inches away from the photocell surface. During the strong visual showers 
would advantageous regroup the photocell units cover the sky area 
the neighborhood 90° from the radiant, where more the radar meteors 
occur. narrowing the apertures both the sensitivity the photocells and 
the accuracy position plotting would proportionally improved. 


The present photocell system, operating under favorable circumstances, 
appears capable detecting first magnitude meteors the average, with 
sensitivity limit about the second third magnitude. The precision time 
measurement the order 0.05 sec. and could improved 
The over-all nonlinear response the system, from the photocells through the 
amplifiers and gates the display C.R.T., precludes the quantitative measure- 
ment the rise and fall the meteor luminosity, though variations lumin- 
osity can detected. 


The value the photocell system aid the visual and radar observa- 
tions has been demonstrated. not unreasonable expect that further 
refinements the instrumentation may put the photocell detector par 
with the human observer respect sensitivity and position plotting, with 
the advantage providing completely automatic permanent record. Direct 
photography will yield more accurate and detailed information about the 
meteor than the photocells can offer, but the time co-ordinate lacking, and 
the sensitivity the existing meteor cameras less than that the present 
photocell system. The Super-Schmidt meteor cameras now under development 
are expected capable photographing third magnitude meteors. Even 
the accurate timing feature only were utilized the photoelectric meteor 
detector might well prove valuable aid the operation the Super-Schmidt 
cameras, because the exact time apparition the meteor required fora 
complete reduction the direct photographs. 
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STUDY IONIC CRYSTALS UNDER ELECTRON 
BOMBARDMENT! 


McLENNAN? 


Abstract 


Electron bombardment experiments have been carried out small crystals 
the alkali halides within the electron microscope. Crystals two size ranges 
were bombarded high intensity, and evidence generalized photographic 
effect within the ionic group solids presented. first group crystal 
specimens ranged size from 0.2 0.002 cm., the bombardment causing the 
formation F-centers and entrapped metal colloids the crystal lattice. 
zation pulses were observed occur the region the specimen during bom- 
bardment. The second group ranged size from 0.01 observations 
the effects bombardment being carried out with electron diffraction tech- 
niques. process has been suggested explain the phenomenon ionic 
crystals which appear lose their centers under high electron beam intensity 
the electron microscope. 


Introduction 


The electron beam employed electron microscopes must relatively 
high energy (order kv. greater) penetrate useful thickness solid 
specimen material. not surprising, then, that such beams high intensity 
(order amp. per sq. cm.) have been observed affect the specimen ma- 
terial. Specimen changes appear based one both the methods 
which electrons may interact with solids: (1) heating, and (2) electronic exci- 
tation. Examples which the former interaction predominates are quite simply 
cited since melting sublimation occurs. the latter case, however, changes 
may caused crystal structure even chemical composition. Watson 
(7) (1948) has discussed specimen changes based variations structure 
iron oxide, tungsten oxide, and tetra-copper calcium oxychloride. Hamm and 
VanNorman (5) (1948) have reported variations structure two metallo- 
organic pigments. The effect temperature these changes discussed 
both authors. Burton, Sennett, and Ellis (1) (1947) have reported inter- 
esting form specimen change the alkali halides. this case clear patches 
are seen form within the crystal after the beam intensity has been increased 
beyond some critical value. The clear patches ultimately spread throughout 
the crystal, leaving behind only thin wall outline its initial position (see 
Fig. 1). Depending upon the alkali halide, one can often observe motion 
material both crystallites and globules inside the crystal during its 
‘evaporation’. Interest this specimen change has prompted the work reported 
this paper. 

Ionic Crystals and Color Centers 


crystals, which the alkali halides are the best examples, are named 
because their atomic components appear exist the crystal lattice ions 
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rather than molecules. important characteristic these materials arises 
from the fact that any crystal defect due missing lattice point (Schottky 
defect) leads holes that bear net charge. These positively and negatively 
charged holes must exist equal numbers electrically neutral crystal and 
give rise the characteristic ionic conductivity and the phenomenon 
color-centers the alkali halides (Mott and Gurney) (6). Defects appear 
exist all crystals, the ratio defects lattice points room temperature 
the best annealed crystals being about 10~*. Color F-centers were investi- 
gated quite extensively Pohl and his coworkers Ger- 
many, the 1920’s. They found that color-centers could produced 
irradiation the crystal with ultraviolet, X-rays, cathode rays, and also 
heating the crystal about 600°C. atmosphere the alkali metal. The 
absorption band giving rise the color crystals the alkali halides 
treated due the transitions electrons trapped missing negative 
lattice points, that is, Schottky defects the positive hole type. the case 
color-centers formed irradiation the crystal with X-rays, Ester- 
mann, Leivo, and Stern (2) (1949) have shown that the density the crystal 
decreases proportion the number color-centers formed. This indicates 
that the defects giving rise color-centers are formed largely the radiation. 
the case crystals heated the metal vapor, Pohl and coworkers have 
shown that excess the metal proportion the number color-centers 
exists the colored crystals. They have also investigated alkali halides colored 
the presence excess metal existing colloid within the crystal lattice 
(Gyulai) (3, 4). this case the color not the same that due color-centers 
nor does appear due electronic transitions. The color seems 
depend, rather, the size the colloidal particles. The color-centers for 
sodium chloride are, for example, yellow, while the color imparted the 
colloid blue-purple. Both forms coloration may produced heating 
the crystal the metal vapor. The color-centers are formed first and are 
the crystal quenched from the high temperature which they 
are formed. the cooling rate decreased the probability forming the 


colloid appears increase rapidly. 


Experimental Investigation 


Electron bombardment experiments were carried out crystals sodium 
chloride and potassium chloride two different ranges size, hereinafter 
called macro and micro. Crystals the size range 0.2 0.002 cm. were chosen 
the macrospecimens since they could observed during bombardment 
under low magnification the electron microscope and, well, examined 
after bombardment under the light microscope The crystals were 
formed from water solution the salt specimen meshes prepared from 100 
mesh per inch copper screen, which was carefully cleaned electrolytic 
polishing. The microcrystals ranged size from and the crystals 
were formed the conventional Formvar-covered nickel stainless steel 
electron microscope specimen meshes vacuum evaporation the salt. 
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Investigation the microspecimens was carried out before and after bombard- 
ment employing electron diffraction techniques. The electron microscope 
was used throughout the bombarding experiments since relatively high inten- 
sity (order 1/3 amp. per sq. cm.) may obtained electron energy 


Bombardment Macrocrystals 


special specimen holder was developed for the University Toronto model 
1944 electron microscope, contain the specimen mesh and one two (as 
desired) insulated electrodes. The electrodes were connected external 
voltage source through high resistance Megohms). Any ionization caused 
the bombardment the region the specimen could then detected 
voltage pulse connecting oscilloscope across the high resistance. 


preliminary electron bombardment experiments, cleaned specimen 
meshes both plain and insulated with thin Formvar film were bombarded 
and attempts were made observe ionization. These blank experiments gave 
indication voltage pulses the oscilloscope; however, indication 
the production secondary electrons was obtained sweeping the fine elec- 
tron beam across one the mesh wires. The steepness the pulse that ap- 
peared depended the sweeping rate. The pulse was obtained with the single 
electrode 250 positive with respect the specimen mesh. Zero voltage 
the electrode showed very slight indication secondary electron col- 
lection while 250 negative showed none whatever. 


With the intense electron beam concentrated crystal sodium chloride, 
pulses positive sign were obtained with negative electrode voltages and pulses 
negative sign with positive electrode voltages. pulses were obtained with 
the collector electrode the same potential the specimen mesh. Calcula- 
tion the ion currents required produce the pulse amplitudes observed 
showed them the order should noted, however, that the 
pulse amplitude was seen depend markedly crystal size, larger crystals 
(0.2 0.02 cm.) giving larger pulses. general induction period from 
sec. was required from the instant playing the concentrated beam 
crystal until pulses were formed. The pulses appeared occur randomly 
and looked like the ionization obtained gas discharge tubes (e.g. 
105, 150, etc.). The same sort phenomenon was observed with potassium 
chloride crystals, but the pulses were general about one-fifth large mag- 
nitude. The use the electron microscope low power facilitated knowing 
whether not the beam was impinging crystal. Pulses were never ob- 
tained with the beam not playing crystal nor were they observed 
certain critical beam intensity was surpassed. was not possible measure 
this critical intensity since the variable was the beam cross section rather than 
the total current. Attempts were made collect the ions produced during 
bombardment, specimen mesh placed the collector electrode. After 
bombardment the collector mesh was examined the methods electron 
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diffraction without removing specimen and collector the air. The results 
were inconclusive since several instances only fine crystallites sodium 
chloride were detected. 


Examination the specimen mesh under the light microscope after electron 
bombardment showed that the crystals sodium chloride were highly colored. 
Crystals which had been irradiated low intensity (considerably less than 0.1 
amp. per sq. cm.) were yellow and amber. However, those upon which the 
intense beam (at 0.1 amp. per sq. cm. greater) had impinged showed con- 
siderable damage and were deep blue-purple color. The latter crystals were 
the ones which also had given ionization pulses under bombardment. Further 
examination the bombarded specimens over period one year has shown 
that the purple coloration remained permanent while the yellow color faded 
within week days. Bombarded specimens potassium chloride showed 
crystals colored both blue and magenta. The former faded within week while 
the latter have remained stable far for vear. 


The unstable yellow coloration bombarded sodium chloride appeared 
due the presence F-centers while the blue-purple color was assumed 
caused the presence sodium colloid within the crystal. check 
this assumption, two simple experiments were carried out purple NaCl 
crystals: (1) purple crystal was dissolved drop water under the light 
microscope, and (2) another purple crystal was sealed damp atmosphere 
chlorine for period six weeks with occasional heating about 200°C. 
the first experiment bubbling was observed until all the blue coloration 
disappeared; the second, examination the crystal after six weeks showed 
some reduction color and the formation new clear crystalline overlayer. 
This was deemed fairly conclusive evidence the presence free sodium. 


Bombardment Microcrystal Specimens 


The Formvar-covered specimen mesh with its layer microcrystals was 
mounted the specially designed specimen holder and placed the electron 
microscope. Bombardment area the specimen defined single mesh 
opening (order 0.010 cm. square) was carried out high intensity after which 
electron diffraction pattern was taken the bombarded material. Speci- 
mens lithium, sodium and potassium chlorides were examined this manner 
both before and after electron bombardment. all cases rings ap- 
peared the diffraction pattern the bombarded salt. further check, 
electron diffraction pictures were taken the same salt with the alkali metal 
evaporated it. The d/n values for the rings produced both 
electron bombardment and the evaporated metal were calculated assuming 
the d/n values the pure salt from X-ray measurements. general the 
lines both patterns may accounted for the presence the 
alkali metal, its oxide, hydroxide carbonate. Tables and are typical 
solution the patterns for bombarded lithium chloride and lithium chloride 
with evaporated lithium metal. The oxides and carbonate appear formed 
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inside the microscope even the relatively high vacuum which exists there 
mm. mercury). Fig. sector (a) shows the diffraction pattern 
lithium chloride itself, sector that caused electron bombardment, and 
sector (c) that due the evaporation metallic lithium the salt 
substrate. Evaporation the metal was carried out inside the electron 
microscope and the diffraction patterns were taken without removing the 
specimen the air. 


The results for sodium and potassium chlorides were similar, the 
rings being caused the presence the free metal, its oxides, hydroxide 
carbonate. well these three salts, four other ionic salts were investigated. 


TABLE 
LiCl BOMBARDED 
Ring No. log d/n d/n log Crystal plane 
14.04 1.1473 2.963 1.6191 LiCt 
14.95 2.78 4438 (1.6185) 2.80 
LiOH: 
16.18 1.2089 2.566 4092 1.6181 LiCl 200 
6 | 2.47 | | Li,COs;: 2. 47 
21.15 1.3253 1.96 2932 (1.6185) 111 
22.89 3598 1.6187 LiCl 220 


TABLE 
LiCl witH EVAPORATED LITHIUM 


Ring No. d/n log d/n log Crystal plane 
15.39 1.1872 2.566 LiCl 200 
110 
21.74 1.3373 1.815 .2589 1.5962 LiCl 220 
211 


log log d/n log K:—based the Bragg law 


where: measured ring diameter, d/n crystal plane spacing, 
the instrument constant. 
4.16 means the crystal plane spacing 4.16 
111 means the 111 plane LiCl. 
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Fic. diffraction patterns (a) normal LiCl, 


(b) bombarded LiCl showing rings, 
(c) LiCl with evaporated indicating some the same rings. 
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These were caesium chloride, sodium bromide, thallium chloride, and calcium 
fluoride, all which showed rings caused the bombardment. Full 
solutions for the rings were obtained the last three cases only, since the 
data for the caesium oxides and carbonates were not available. Again the 
presence the metal its oxidation products could account satisfactorily for 
the extra rings. 

Discussion 


appears from this work that something comparable the photographic 
process occurring all the ionic materials investigated, the 
the silver halides appearing part more general phenomenon 
characteristic ionic crystals. For this reason the process may termed 
“generalized photographic 


The following process now suggested explanation the specimen 
change observed occur the alkali halides under electron bombardment. 
First all the formation color-centers appears equilibrium process. 
This was observed Estermann, Leivo, and Stern (2) (1949) when they 
pointed out that X-ray irradiation for periods longer than hr. failed 
show any increase the number color-centers. may then write 
equilibrium equation describe the formation color-center: 


F-center 
(A) (B) (C) 


(A) represents ionic bond between sodium and chlorine ions the sodium 
chloride lattice; the presence radiation (in this case, the electron beam) 
will ionize one the lattice points. The excited lattice point may then capture 
the most loosely bound electron the Na-Cl configuration. This the valence 
electron sodium which has moved into the valence band form 
and which accounts for the ionic binding. (B) intermediate configura- 
tion which sodium ion and neutral chlorine atom are left with very little 
binding between them. The chlorine atom then assumed move off inter- 
stitially under the temperature vibrations the lattice. The sodium ions 
surrounding vacant negative ion lattice point (Schottky defect) cause 
positive hole which may then trap free electron form F-center (C). 
The presence free chlorine the lattice may cause this process the 
other way and re-form normal lattice. the intensity the bombarding 
beam increased, the equilibrium moves the right; any localized increase 
the amount free chlorine will cause the equilibrium move the left 
again, thus reducing the number F-centers. the number F-centers 
increased increasing the lattice will become more and more unstable 
until collapses. However, increase the amount free chlorine and 
the temperature the lattice from electron bombardment will cause the 
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chlorine escape the outside. There thus tendency have more free 
chlorine the outside the crystal than the center; logical then 
that when collapse does begin should begin the center. This catastrophic 
degradation the crystal proposed the reason for the specimen change 
observed Burton, Sennett, and Ellis (1), for the blue-purple macrocrystals 
and for the rings the diffraction patterns bombarded microcry- 
stalline specimens. 

Since the catastropic degradation seems occur only after critical inten- 
sity has been surpassed, would interesting whether intense beam 
soft X-rays could made degrade the alkali halides similar fashion. 
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ENERGY SPECTRUM NEUTRONS 
FROM 


Abstract 


The energy spectrum the neutrons emitted when beryllium bombarded 
with alpha particles from C’) has been determined measuring the 
ranges and directions recoil protons produced cloud chamber the 
neutrons. The maximum energy the neutrons was found about 15.5 
Mev., within experimental error the calculated maximum, 14.4 Mev. Many 
neutrons were found with energy about 2.15 Mev. and the structure the 
remainder the spectrum can attributed the existence energy levels 
about 7.1 and 4.5 Mev. angular distribution the protons scattered 
the neutrons from this source was isotropic within experimental error. 


Introduction 


data have been found the energy distribution the neutrons from 
thorium C’) beryllium source. have determined this distribution 
measuring, cloud chamber, the energy transferred the proton the 
neutron-proton scattering interaction. 


have investigated the angular distribution the recoil protons and have 
found isotropic for the neutron energies found here. This agreement 
with the distributions found many other investigators (1, 8). 


Apparatus 


The cloud chamber used was the rubber diaphragm, 
type, in. diameter and in. depth. This chamber has been described 
detail Lillie (6). Stereoscopic photographs were taken with two Kodak 35, 
mm. cameras rigidly mounted above the chamber aluminum casting. 
Illumination was provided discharging 100 condensor charged 2000 
v., through two G.E. FT-126 mercury vapor discharge lamps mounted 
opposite sides the chamber. 


The operation the complete cloud chamber unit was controlled series 
thyratron delay circuits cascade and was completely automatic. This 
control scheme was chosen preference the motor-driven cam system 
because its greater flexibility. 

The reprojection unit used studying the track photographs 
described elsewhere (2). The track images were reprojected through their re- 
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spective camera lenses. When the two images track were brought into 
coincidence screen all distortion due asymmetry the cameras their 
mounting was automatically avoided. 

The neutron source used shown Fig. This source was designed that 
could placed entirely within the cloud chamber. The aluminum disk 
carrying the active deposit mm. diameter and the beryllium 
has areal density 0.932 gm. per 


Fic. shows the neutron source with the glass walls the cloud chamber the rubber 
gasket for sealing the chamber the brass cylindrical sleeve the mm. lead sheet the beryllium 
metal sheet the aluminum disk carrying the thorium active deposit and the removable plug 


which holds the source position 


The mm. lead sheet served reduce the beta and gamma background 
the chamber due source activity. This lead affected the energy less 
than the neutrons passing through with energies greater than about 
Mev. since the mean value for the total cross section these neutrons 
lead about 5.5 barns (5) corresponding mean free path about cm. 
The source was mounted inside the glass chamber wall shown and was 
constructed that the aluminum disk could removed for reactivation without 
disturbing the gas mixture the chamber. This construction also prevented 
the gas mixture from coming into contact with the thorium C’), thus 
preventing contamination the chamber. 


Experimental Method 


The thorium preparation used obtaining the thorium active deposit 
sources was intimate mixture radiothorium and ferric oxide carrier. This 
mixture was placed brass cylinder with insulating cover such that the 
alumimum disk could suspended the space above the radiothorium and 
maintained 400 negative potential relative it. Strong alpha-active 
sources were obtained this way with even half-hour exposures. com- 
parison with standard radium source the activity typical deposit has 
been estimated millicuries. The sources were removed from the charging 
device least one-half hour before use. Hence the alpha particles were those 
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was not possible operate the chamber which was available for this work 
pressures greater than atmospheric. Since the stopping power hydrogen 
nitrogen and oxygen atmospheric pressure the noncondensable gas order 
stop higher energy particles within the chamber and make the conditions 
for good tracks less critical. This mixture simulates standard air. 


Fic. This shows typical proton recoil. The neutron source shown position the 
chamber wall the lower right. The the serial number assigned to, and photographed with, 
this exposure. 


This gas mixture was then saturated with 50% ethyl alcohol, 25% water, 
and 25% acetone mixture, with the result that the number hydrogen atoms 
available the chamber was between and 60% the number when using 
pure hydrogen. 


Over 1000 stereoscopic exposures were taken Ansco Ultra Speed Pan film 
and were given standard processing Kodak developer and Edwal 
Rapid Fixer. Fig. shows typical proton track. 
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Before and after each series photographs polonium source was inserted 
the chamber and the alpha tracks photographed for use reducing the 
observed proton ranges ranges standard air. assume that the stop- 
ping power the gas mixture relative air constant all points along the 
range and the apparent range alpha particle measured under con- 
ditions identical with those during the measurement the apparent proton 
range then may write LaX where the known range 
the alpha particle standard air and the reduced value referred 
standard air. 


using the reprojection system make measurements the track photo- 
graphs (2) the track images were projected into the object space corresponding 
the chamber volume and onto plane, white, circular screen. This screen was 
magnetically clamped such way that was free rotate about any dia- 
meter passing through its center while this center could moved any point 
the object space. making measurements the proton tracks, the track 
images were brought into coincidence the screen surface; the screen was 
moved bodily until this track lay diameter and the screen support was 
locked rigidly the base plate the unit energizing the base magnet. The 
screen proper was then rotated about this diameter until the images the 
neutron source were brought into coincidence and then was locked place 
energizing the support magnet. 


Proton ranges were measured with mean error less than mm. and the 
angles between the track and line joining the beginning the track the 
center the neutron source were measured the nearest degree. Since this 
source had finite dimensions was necessary reject all tracks starting within 
the source order reduce the probable errors angle. From the 
geometry the system, using average distance cm. from the source 
the beginning the track, the mean expected error the angle recoil 
less than two degrees. All tracks which struck the chamber walls were rejected 
the energy measurements but not the measurement angle the clear 
length was greater than cm. 


When the angles and ranges tracks indicated neutron energy far excess 
any expected from the source (25 Mev.) the tracks were rejected. These 
were doubtless due neutrons from the source which had previously been 
scattered the material the chamber. The number such tracks was less 
than the total number measured. 


The measured proton ranges were converted the reduced range referred 
standard air making measurements the ranges polonium alpha 
tracks photographed before and after each series proton track exposures 
and using the conversion relation noted previously; the corresponding energies 
were obtained from the curves given Livingston and 
Bethe (7). 


Conservation energy and momentum leads the well known expressions 
E,= and tan where the proton recoil energy, 
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the incident neutron energy and the angle recoil measured between the 
proton track and the initial direction the neutron. Using this relation the 
neutron energies corresponding each proton track were calculated. The 
energy distribution obtained this way must corrected for the variation 
the scattering cross section with energy order obtain 
the true energy spectrum. These corrections were obtained from the smooth 
curve given Goldsmith, Ibser, and Feld (5), which gives the measured 
variation the total scattering cross section with energy. Fig. the final 
corrected spectrum based measurements 750 recoil protons. 


a 


So 


NUMBER NEUTRONS (Arb. Units) 


NEUTRON ENERGY 


intervals. 


The angular data obtained the energy measurements were used check 
the angular distribution the recoil protons. Since this distribution has been 
studied many investigators was considered interest compare our 
results with those obtained previously. 


Let represent the element solid angle and the angle recoil dw; and 

the center mass system co-ordinates and and the labor- 

atory system. the center mass system the number protons recoiling 

through unit solid angle angle recoil may represented 
function this angle, 

dw, 


From simple trigonometry and 
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definition the angular distribution isotropic constant 
against should have flat maximum 45° and should drop zero 


LEGEND: ANGLE RECOIL 


BROMLEY CHAMPION and POWELL 


Fic. experimental isotropic distribution. results are quoted 
Rosenfeld (8).) 


The ratio the observed number tracks each five degree interval 
the number calculated for isotropic.distribution was computed and plotted 
against the mean angle the interval. The result shown Fig. together 
with the results four other investigations, (1), (3), (4), (8), whose data have 
been reduced this form for comparative 


Discussion Results 


noted previously, the only alpha particles are those from C’) 
6.0537 and 8.7759 Mev. neglect the very rare long range alphas from the 
thorium C’. 


These alpha particles penetrate varying distances into the beryllium crystal 
structure before initiating neutron producing reaction. Furthermore, the 
angle neutron emission unknown and the amount momentum carried 
off the product nuclei varies over wide range. For these reasons should 
expect find, especially source which the alpha emitter and beryllium 
are intimately mixed, that the neutron energy groups would combine give 


= 
uJ + b ee 
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continuous spectrum. With the source use here the incident alphas are 
more nearly monodirectional impact, whence should expect that the 
more intense groups, this case, combinations two groups, might appear 
superimposed continuous spectrum. The spectrum obtained (Fig. shows 
these features. The only maxima which are statistically justified con- 
sidering are those corresponding neutron energies about 14.5, 11.5, 7.3, 
and 2.15 Mev. 


The major neutron producing reaction this source (a, C}?. Tol- 
lestrup, Fowler, and Lauritsen (9) give 5.708 Mev. the value. These 
authors also give value —1.67 Mev. for the reactions (y, 
which are also energetically possible this 
source. The cross sections for these reactions are low and this type source 
the reaction particular unimportant since the beryllium which 
the neutrons are produced the thin sheet. 


The calculated maximum energy for source this type 14.49 Mev. for 
the 8.78 Mev. alphas and 11.76 Mev. for the 6.05 Mev. alphas. 


The angular distribution recoil protons indicated 45° the most probable 
angle recoil and have previously given the mean expected error the 
angle less than two degrees. This gives 0.069 E,, where the 
neutron energy. 15.5 Mev., dE,= 1.1 Mev., hence the observed 
maximum energy about 15.5 Mev. within the experimental error the 
calculated value. 


not consider that the statistics are sufficiently good warrant cal- 
culation the neutron distributions for transitions the excited states 
was done detail Whitmore and Baker (10). note, however, that 
the maxima 11.5 and 14.5 can explained considering that they are 
due transitions the ground state which the full alpha energy 
has been transferred the neutron. similar reasoning can consider that 
the maximum 7.3 Mev. due the superposition two groups, one cor- 
responding transitions the 7.1 Mev. excited level and resulting from 
the 8.78 Mev. alphas and the other corresponding transitions the level 
about 4.5 Mev. resulting from the 6.05 Mev. alphas since have 14.49 
Mev. and 11.76 Mev. 


The relatively large number high energy neutrons found here considered 
due the particular type source used which the mean angle neutron 
emission with respect the direction the incident alphas much smaller 
than source consisting mixture the alpha emitter and beryllium. 


The large number neutrons the 2-3 Mev. range can considered 
part the previously mentioned continuous distribution. The sharp decrease 
below Mev. considered due experimental cutoff. For neutron energies 
below Mev. only the recoil protons making small angles with the neutron 
momentum would have sufficient energy detected owing the relatively 
high background the chamber. 


™ 
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The statistics here are not sufficient enable give any definite evidence 
the existence nonexistence the previously found energy levels the 
nucleus about 7.1, 4.5, and 2.5 Mev. Whitmore and Baker (10) have 
summarized the latest available data these levels. 


seen from Fig. the angular distribution the recoil protons isotropic 
within the precision this work. This precision not great enough indi- 
cate any small deviations from isotropy, they exist. The excess 
the number tracks low angles can explained noting that for these 
small angles recoil greater fraction the neutron energy transferred 
the protons. The probability detection thus increased. For large angles 
the energy transferred small and the probability detection above the 
chamber background decreased, thus explaining the experimental deficit 
the number tracks large angles. Many investigators have studied this 
distribution and have found conclusive evidence for any departure from 
isotropy the energy range considered here. 
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PULSER FOR CYCLOTRON 


Abstract 


pulser for the R.F. oscillator the McGill synchrocyclotron described. 
This pulser triggered selected precise frequency and controls the oscillator 
for selected pulses for predetermined exact periods continuous operation. 


Introduction 


the direct exposure photographic plates studies very short 
activities frequently necessary operate the cyclotron for short periods. 
order that such periods reproduced with accuracy pulser has been built 
for controlling the R.F. oscillator the McGill synchrocyclotron. Previous 
pulsers reported have either used mechanical device produce the triggering 
pulse (3), have shut off the oscillator for certain period each f-m cycle 
grid blocking (1). the operation the McGill pulser the oscillator 
normally not oscillating, and accurately triggered the desired frequency 
the f-m cycle the absorption power from low power oscillator the 
dead tank circuit. 

Description Apparatus 


The frequency modulated voltage for the McGill cyclotron dee provided 
class grounded grid oscillator employing 880 triode (2). Since the 
plate the oscillator tube requires continuous water cooling, the tube con- 
nected with the plate d-c. ground and the filaments high negative poten- 
tial. class operation the grid the tube is, the average, several hundred 
volts negative with respect the cathode. The d-c. supply capable 
delivering kv. and thus the grid may low kv. below ground. 


The operation the oscillator such that the actual duty cycle comprises 
only one-sixth the total f-m cycle. view this was considered advisable 
arrange the pulser such way that the oscillator normally not oscil- 
lating, and permitted oscillate only during the period the duty cycle. 


With this end view, tests were made determine the best method 
blocking the oscillator means sufficiently large negative voltage applied 
the grid, the same time keeping mind that any switch inserted this 
grid circuit must present the proper impedance the tube when operating 
over the duty cycle. was found that two 304TH triodes, connected parallel 
and series with 500 resistor, fulfilled the latter condition. kv. power 
supply was inserted across the 304TH combination with connections made 
such way that when the are nonconducting the grid the oscil- 
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lator tube biased negative this amount. Under these conditions the 880 
tube will not oscillate. The 304TH combination then forms effective switch 
since they can made nonconducting and switched during the duty 
cycle. This circuit shown Fig. 


R.F. OSCILLATOR 


R; = 13k Ry = 10k C3 = 0.1 mf. 
R, = 22 k Rio = S7 k 

R; = 2 Meg. Ru = 390 k 


The are biased cutoff current flowing through the plate re- 
sistor the 6L6. The voltage for this tube and its control unit provided 
300 power supply whose positive side connected the negative 
terminal the kv. supply. The unbiased grid the 6L6 capacitance- 
coupled the plate the normally nonconducting tube uni-vibrator. 
Thus, when the uni-vibrator its stable state, the grid the 6L6 the 
same potential the cathode and hence the tube conducts strongly, providing 
sufficient bias its load resistor keep the cutoff. However, when 
the uni-vibrator triggered into its unstable state, large negative square 
pulse appears the grid the 6L6 and ceases conduct, causing its plate 
rise the potential the 300 supply, and hence removing the grid bias 
the Under these conditions the are capable conducting 
and the R.F. oscillator proceeds oscillate. The duration the oscillation 
controlled the length the square pulse appearing the grid the 6L6, 
which turn governed the circuit constants the uni-vibrator. These 
may adjusted give any length pulse desired. 


04TH 304TH 880 
O+ 
2000V 
PULSE 
TRANS. 
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The grid the uni-vibrator potential approximately 11,000 
below ground, and since the triggering pulse must come from the ground level, 
the most suitable arrangement pulse transformer capable withstanding 
this voltage. 


The duty cycle the oscillator extends from 25.2 Mc. 21.9 Mc. whereas 
the frequency modulation cycle from 27.9 20.9 Mc. From the point view 
other circuits, such the timing pulser for the ion source, which must 
pulsed the start the duty cycle, desirable have oscillator 
operating slightly before 25.2 Mc. 


Ry = 100 k 


order obtain triggering pulse for the pulser the circuit shown Fig. 
employed. local low power oscillator coupled with tuned circuit 
acting detector, which coupled turn the main tank circuit. Both 
the local oscillator and detector are tuned Mc. Each time the main tank 
sweeps through this frequency, absorbs energy and the amount absorbed 
depends upon the coupling. The drop level the detector fed into 
cathode follower which feeds the pulse the delay and gating unit shown 


Fig. 


The pulse obtained from the detector occurs 27.0 Mc. and necessary 
delay this pulse until the frequency the main tank circuit reaches 25.4 Mc. 
This pulse, amplified the 6SJ7, triggers the 2050 thyratron, whose output 
triggers uni-vibrator and the square pulse obtained differentiated give 
positive pulse, which fires the second 2050 thyratron. The pulse from 
this thyratron can fed directly into the pulse transformer, thus triggering 
the main pulser. The length the square output pulse from the uni-vibrator 
adjusted that the delay the appearance the positive pulse sufficient 
allow the main tank circuit reach 25.4 Mc. Since this time depends 
directly the speed the rotary condenser, the delay must variable over 
fairly wide limits accommodate the speeds which may employed. 


00 
=) =) 
c 
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considerable amount energy absorbed the detector when the main 
oscillator on. prevent this from affecting the control circuits the 6SJ7 
gated. The output pulse from the first 2050 triggers second uni-vibrator. 
The negative output pulse from this uni-vibrator applied the grid 
normally conducting 6SJ7. The voltage drop the cathode resistor this 
tube normally provides the screen voltage for the 6SJ7 amplifier. However, 
when the cut off, the screen the 6SJ7 amplifier drops zero and 
the tube ceases conduct. this manner any disturbance attempting cause 


spurious triggers for the main oscillator blocked. 
Rig 


INPUT 


R; = 100 k Ris = 2k C3 => 1 mf. 
R, = 100 k Rig = 22 k Cs = 0.001 mf. 
Run = 47 k Rog = 25 Meg. 
Rus = 15k (10 w) Rog = 15 Meg 
Ris = 47k Rx = 1k 


Owing slight unevenness the teeth the rotary condenser, not 
possible tune the local oscillator the exact top end the modulation 
cycle. Hence the main oscillator tank circuit sweeps through the frequency 
the local oscillator twice each cycle, two pulses appearing the input the 
6SJ7 amplifier. The first pulse passes through and triggers the first 2050 thy- 
ratron, and this turn triggers the gate which blocks the second pulse, well 
‘hash’ from the main oscillator. 


OO 
Ris 
q 
+300 
R, 
R 
R, 
20 Ros 
0 TO -150 
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Since has been found necessary some cases (i.e., nuclear plates) expose 
for very short periods, the circuit shown Fig. can inserted between the 
delay unit and the pulse transformer. The output pulse from the delay unit 
thyratron applied the grid second thyratron, which shares its cathode 
resistor with 6F6. The bias this thyratron adjusted that, with the 
6F6 conducting hard providing additional bias the cathode the 2050, the 
tube does not fire pulses from the delay unit. However, when the 6F6 
made nonconducting this additional bias removed and the thyratron fires 
the desired pulses. This condition achieved allowing condenser pre- 
viously charged 300 discharge itself through the grid resistor the 
6F6, thus cutting off. The duration the cutoff determined the time 
constant this condenser and grid resistor. The latter variable provide 
for passage pulses for times from zero one second. 


+300 


The flexibility the control units such that any manner triggering may 
accomplished inserting suitable circuits place the circuit Fig. 
addition this, all the controls for these circuits can operated the 


control desk the cyclotron. 
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THE CALIBRATION PLATINUM RESISTANCE 
THERMOMETERS THE TEMPERATURE RANGE 11° 


Abstract 


set six platinum resistance thermometers form suitable for low 
temperature calorimetry has been made and calibrated the region 11° 90°K. 
intercomparison with similar thermometer which had been calibrated 
the National Bureau Standards. Above 90°K. calibration has been made 
the International Temperature Scale. 

Using the intercomparison data, has been possible derive method 
whereby for these thermometers the scale for the region 20° 90°K. may 
found within 0.002°C. means fixed points. The method applies 
function’ the type used the National Bureau Standards (13), plus 
corrective term which depends upon the resistance the thermometer the 
boiling point hydrogen and upon the normal constants which are determined 
for the International Temperature Scale above 90°K. 


Introduction 


preliminary step the construction low temperature calorimeter 
assembly, was necessary make and calibrate set platinum resis- 
tance thermometers. view the fact that the measurement the temper- 
ature has important bearing the precision claimed for the calorimetry, 
seemed desirable describe the thermometers and manner which they 
were calibrated. addition, factors have been found which pertain the 
general problem the establishment temperature scale the region below 
the boiling point oxygen. 


The problem establishing temperature scale for the region below 90°K. 
course well recognized (4). the National Bureau Standards pro- 
visional scale has been established (7) terms gas thermometer, with 
view seeking first all uniformity temperature scales for the various 
laboratories doing low temperature work. That this much might accom- 
plished, series thermometers has been made the Leeds and Northrup 
Company, and calibrated the National Bureau Standards. One these, 
#718154, was purchased the National Research Council and serves the 
standard for the calibration working thermometers used this laboratory. 


The objective further work this temperature scale the establishment 
scale which agrees closely with the thermodynamic scale and which may 
realized terms the resistance platinum thermometers through the 
use fixed points and suitable interpolation formula. scale established 
this way could serve extension the present International Temperature 
Scale (12). Perhaps the greatest difficulty has been obtain suitable inter- 
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polation formula. empirical function’ the type first proposed 


Roo. 190 — Rss. 363 


has been tested out the National Bureau Standards. The same tabulated 
relation between and used for all thermometers. Such function may 
used for the region between the triple point oxygen (54.363°K.) and the 
boiling point oxygen (90.190°K.) with considerable success, but similar 
Roo. 273 

where the boiling point equilibrium hydrogen taken 20.273°K., may 
give deviations large 0.05° from the real temperature. More recently, 
Hoge (6) has suggested another method which linear combination two 
different functions used. the basis the present work, has been 
possible derive method which uses single function’ and corrective 
term which dependent only upon knowledge the resistance the un- 
calibrated thermometer certain fixed points. While the method operates 
very satisfactorily the present instance, the thermometers here described 
have characteristics which are not too different, hence would seem desirable 
for the method tried with thermometers made other laboratories. 


Experimental 


Description the Thermometers 


The six resistance thermometers were made following the method Meyers 
(9), the two designated T-4 and T-5 being made first some nine months 
before the other four, which were taken from set made the one time. 
All the thermometers were made from the same batch wire obtained from 
the Sigmund Cohn Company New York. The dimensions the mica crosses 
were made such that the thermometers fitted closely inside platinum thimbles 
mm. long, mm. diameter, and 0.2 mm. wall thickness. After being annealed 
constant ice point resistance heating temperatures above 450°C. 
air, and after the steam and sulphur points had been determined, the thermo- 
meters were sealed off with soft glass cap, the four leads (two current and two 
potential) being sealed through the cap. Before sealing, the thermometers were 
filled with amount dry helium corresponding about 1/3 atm. pressure 
room temperature. 

The Leeds and Northrup thermometer 718154, hereafter designated 
T-S, had been made similar way but somewhat smaller dimensions, the 
‘platinum thimble being approximately mm. long, and mm. diameter. 


Determination the Fixed Points the International Temperature Scale 


For the determination the fixed points the boiling points water and 
sulphur, hypsometers the type described Beattie (1, were used. The 


ail 
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atmospheric pressure was determined with mercury manometer conven- 
tional design, the heights the mercury columns being compared with 
calibrated steel scale mounted the arms the manometer. Appropriate 
corrections were made for capillary depression the mercury, and the pres- 
sures were reduced the beight mercury column and standard 
gravity the level the resistance thermometer coil. From large number 
measurements, appeared that the accuracy pressure determination was 
0.03 0.04 mm. 


The ice point and oxygen point determinations were made the conven- 
tional manner, the oxygen point apparatus being essentially that described 
Scott (11). 


Resistances were determined with White double potentiometer. For the 
fixed point determinations, steady current about 1.2 ma. was passed 
through the thermometer circuit. The sensitivity the measuring system was 
such that change 0.00001 ohm could detected. For most the ther- 
mometer intercomparisons which are described below, steady current ma. 
was used. The current was raised ma. for some the intercomparisons 
below 30°K. 

TABLE 


INTERNATIONAL TEMPERATURE SCALE CONSTANTS 
FOR THE RESISTANCE THERMOMETERS 


| 
Average Average Average 


Unrounded values determined the National Bureau Standards. 
the basis later intercomparisons adjusted the equivalent 0.008° this value. 


Table are given the values the constants the thermometers the 
Dusen equation 


1007 \100 100 100 


columns and are given the average deviations the determinations 
and respectively. The numbers brackets indicate the number 
individual measurements which were averaged. 

Intercomparison the Resistance Thermometers the Region below 90°K. 


For the intercomparison, the thermometers were set into re-entrant wells 
copper block with Woods metal. Four wells were provided the thermo- 


aid 
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meters were compared with T-S three time. The copper block was hung 
within the adiabatic shield low temperature calorimeter assembly which 
has been described elsewhere (10). 

Altogether, some comparisons each thermometer with were made 
the temperature range 11° 90°K. keeping the difference temperature 
between the block and the adiabatic shield small possible, practically all 
the intercomparisons were performed without the temperature the block 
changing. any one set intercomparisons, the measurement the resis- 
tance was usually made alternately with the measurement the resis- 
tance each the other three thermometers. 

The adiabatic shield had been provided with vapor pressure thermometer 
bulb, and this was used check the indication the boiling point 
equilibrium hydrogen. The bulb was partially filled with about gm. 
which was prepared the manner described Taylor Diamond (14) 
and which served the catalyst for converting the hydrogen the equilibrium 
mixture. About 0.04 mole pure hydrogen was condensed into the bulb and 
after period about three hours the steady value the pressure and the 
resistance were determined simultaneously. The indication the 
sensitive thermocouple between the block and the adiabatic shield was also 
recorded. The temperature the vapor pressure thermometer was calculated 
from the formula Woolley al. (15) for equilibrium 


was found that the indications the thermometer and the vapor pres- 
sure thermometer agreed 0.008°. The vapor pressure thermometer bulb 
being situated the adiabatic shield was really not ideal for the experiment, 
the agreement considered very satisfactory. 

From the intercomparison data calibration table for each thermometer, 
T-1 was obtained use method suggested Hoge (5). For the 
smoothing the data the resistance uncalibrated thermometer was ex- 
difference AR’ plotted against with AR’ being made 
zero about 26° and 86°K. The scatter the points plotted from the 
smooth curves was general less than 0.00005 ohm. the same graph was 


temperature table for each thermometer directly, with intervals such 
permit linear interpolation. For each value preliminary value the 


temperature was obtained through and the calibration table 


T-S. The correction, get the final temperature was determined from 


AR’ 


= 
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The scale derived was joined smoothly with the International Temperature 
Scale the oxygen boiling point adjusting the value for each thermo- 
meter that the value AR’ Reatc. 90.19°K. fell the 
smooth curve the AR’ against plot. Such adjustment permissible 
because the intercomparisons may performed more precisely than the deter- 
mination the oxygen point for the International Temperature Scale. For five 
the thermometers the value was changed the equivalent 0.001° 
less. For adjustment equivalent 0.008° was required, indicating 
that systematic error had been made the oxygen point determination for 
that thermometer. However, view the consistency the behavior the 
other thermometers, did not seem worthwhile redetermine the oxygen 
fixed point for 

While the scales are joined smoothly 90.19°K., continuity the deri- 
vative resistance with temperature does not follow necessarily. The deri- 
vative can discontinuous owing deviation the International Tempera- 
ture Scale from the thermodynamic scale above the boiling point oxygen. 


Method Determining the Scale for the 
Region 20° 90° means Fixed Points 
R Rss. 363 


The proposal Cragoe (13) use the function 
Roo. 190 Rss. 363 


was tried with the thermometers used the present work.* 


The resistances the thermometers the oxygen triple point (54.363°K.) 
were taken from the intercomparison data and not determined directly. Table 
shows the maximum error arising from the determination the temperature 
using this function rather than the intercomparison data. 


TABLE 


MAXIMUM ERROR INVOLVED USING FUNCTION 
THE RANGE 54.363° 90.190°K. 


Thermometer max. degrees 


T-S 0.0002 
T-1 
T-2 
T-3 
0.0011 
T-5 0.0010 
T-6 


The small value the errors shown the table suggests that these thermo- 
meters have similar characteristics the one which the table was based. 


attempt was now made determine the temperature scale for the region 
from the boiling point equilibrium hydrogen (20.273°K.) the boiling point 


copy the appropriate table relating the function the temperature was placed our 
disposal the National Bureau Standards. 
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oxygen (90.190°K.) using single function. new table was calculated 


(Table with the function based upon T-S. The 


Ro0.190 Roo. 273 
value was taken from the National Bureau Standards calibration. 
For the thermometers plot was made against for three 
drawn through the points. These curves are shown solid lines Fig. where 


~.0002 


~.0002 


~.0002 


20 30 40 60 7° 60 90 


50 


Fic. relationship between and temperature for the six thermometers. 


determined from the intercomparison data. 
from Equation (IV). 
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seen that for and T-5 deviation large 0.009° found about 

means was now sought obtain general formula express the observed 
discrepancies terms the thermometer constants shown Table Empi- 


TABLE III 
Z-TEMPERATURE RELATIONSHIP FOR THE FUNCTION 
Ro0.190 Roo.273 

|| em 
—0. 003680 409139 


rically was found that the curves could approximated most simply 


The quantities and were determined from the following conditions 
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(c) = = 0.01816486 (Roo. 190 Roo. 273) 


where Ro, and are the constants the Callender-VanDusen equation. 
Equation derived the assumption that 90.190°K. equal 


the slope obtainable from Rinternat. 


The factor 0.01816486 was found from the requirement that 


90.190°K. 


From the above, Equation (I) becomes 


sin (2.5745T 52.19), 


where the phase 52.19) expressed angular degrees. Using 
Equation (III) the solid curves shown Fig. could approximated quite 
reasonably, the largest deviation being one 0.003° for thermometer 


Closer examination the data showed that except for the slopes 
the versus curves 90.190°K. departed significantly from those deri- 
vable from the equation. Blaisdell and Kaye (3) have 
pointed out that part the internal inconsistency the International Tem- 
perature Scale can due differences for various thermometers. Further, 
from work Hoge and Brickwedde (8) appears that the indications 
platinum resistance thermometers are generally too low between 
and 0°C. for large values One might expect then that the observed differ- 
ences slope the oxygen point would dependent upon was found 
that the differences could reduced considerably replacing Equation 


sin (2.5745T 52.19). (IV) 


The plots versus for each thermometer determined from 
are shown dotted curves Fig. The deviations the calculated curves 
from those derived from the intercomparison data are nowhere larger than 
ohm, terms temperature deviation nowhere larger than 0.002°. 


summarize, proposed that the use the table (Table with 
appropriate corrections evaluated Equation (IV), may provide con- 
venient means accurately establishing the temperature scale for platinum 


— 
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resistance thermometers the region 20° 90°K. The addition but one 
fixed point (at the boiling point hydrogen) those now required for the 
International Temperature Scale necessary. The generality the method 
can only evaluated after has been tested with thermometers character- 
istics somewhat different from those used this investigation. suggested 
that the scale the region below 20°K. may established most conveniently 
direct comparison the platinum resistance thermometer with hydrogen 
vapor pressure thermometer. The vapor pressure thermometer would re- 
quired any case for the determination the fixed point the boiling point 
hydrogen. 
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THE RAMAN SPECTRUM FLUORINE! 


ANDRYCHUK 


Abstract 


The rotational and vibrational Raman spectra gaseous fluorine atm. 
pressure were obtained. From these spectra the constants the molecule 
the electronic ground state were determined be: 

0.8828 
The intensity alternation the rotational spectrum shows that the nuclear spin 


fluorine and that the nuclei follow Fermi statistics. 


Introduction 


Spectroscopic data for the determination the molecular constants 
fluorine the electronic ground state are incomplete. The molecule does not 
exhibit infrared absorption spectrum since has dipole moment (9). 
The emission spectrum observed Gale and Monk (6) well Aars (1) 
electric discharge through fluorine yields data about the electronically 
excited states 'II and For these states the internuclear distance was found 
1.282 and 1.475 respectively, and the corresponding vibrational 
frequency 1139.8 and 977.4 Wartenberg, Sprenger, and 
Taylor (15) showed that the absorption spectrum fluorine the visible and 
near ultraviolet regions, unlike the absorption spectra chlorine, bromine, 
and iodine, entirely continuous and does not therefore lead any infor- 
mation about the vibrational and rotational constants fluorine the 
ground state. 

approximate value 1.435 for the internuclear distance the ground 
state fluorine was obtained Rogers, Schomaker, and Stevenson (14) 
means the electron diffraction method. Using earlier electron diffraction 
value 1.46 obtained Brockway (4) and applying rule (3) 
Murphy and Vance (12) estimated the vibrational frequency fluorine 

The only direct spectroscopic method determining these quantities would 
appear either study the Raman spectrum study the absorption 
emission spectrum the far ultraviolet. Neither these methods has been 
successfully applied until now. unsuccessful attempt observe the Raman 
spectrum has been described Garner and Yost (7). Considerable improve- 
ments light sources and the speeds photographic emulsions made 
seem worth while make renewed attempt observe the Raman spectrum 
gaseous fluorine. already mentioned preliminary note (2) the vibra- 
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tional Raman spectrum fluorine was indeed found and recently the rota- 
tional spectrum has been obtained and resolved. The present paper gives 
more detailed report this work. 


Experimental Procedure 


The Raman tube was patterned after the design developed Crawford and 
Welsh (5) the University Toronto. Fused quartz tubing was employed, 
since according Leech (10) quartz especially resistant dry fluorine free 
from hydrogen fluoride. The Raman tube was in. diameter and in. long, 
and consisted three sections shown Fig. The front section AB, in. 


Fic. diagram the apparatus. 


long and painted black the outside contained copper diaphragms which shield- 
the window the Raman tube from direct light. These diaphragms were 
matched optically the spectrograph means the lens The center sec- 
tion was the irradiated portion and was in. long. The third section CD, 
in. long and painted black the outside, was the horn end the Raman 
tube. cell containing dust free sodium fluoride pellets was fused the horn 
end the tube. 


The tube was baked under vacuum remove moisture from the walls and 
then sealed off from the vacuum systems three liter Florence flask con- 
taining the dry fluorine was attached the Raman tube means the ground 
glass joint The fluorine was stated 97% pure. The space was 
evacuated first, then sealed off order fill the Raman tube with 
fluorine, the break-off seal was broken and the two stopcocks were opened. 
Any hydrogen fluoride which may have been present was removed from the 
fluorine passed over the sodium fluoride pellets When the tube was 
filled the stopcocks were closed. 


The ground glass joint and the stopcocks were made air tight with the aid 
fluoro-carbon grease. This method handling the fluorine the Raman tube 
prevented the clouding the tube for long one week. 


| 
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The light source consisted two water-cooled mercury arc lamps the 
type developed Crawford and Welsh the University Toronto (5). 
description these lamps given Rank and Kagarise (13). The lamps were 
operated with power input kw. Magnesium oxide coated reflectors were 
used with the lamps, since Menzies and Skinner (11) have shown that the use 
this type reflector greatly increased the light intensity. 


Zeiss three prism spectrograph was employed. achromatic doublet lens 
having focal length 50.8 cm. and aperture F/10 was adapted the 
camera lens. The linear dispersion was 4047 The rotation 
and vibration spectra were photographed Eastman Kodak plates. 
Two plates were obtained using spectral slit width 1.2 and ex- 
posure hr. third plate was obtained using spectral slit width 0.8 


The spectrum iron arc was photographed juxtaposition with each 
Raman spectrum. The positions the Raman lines and the calibration lines 
were read Hilger comparator one thousandth millimeter. 


For each plate the deviations the comparison lines from linear rela- 
tionship between line position and wave number were plotted and the wave 
numbers the rotational lines were determined with the aid this correction 
curve. The values for the displacements from the exciting line given Table 
are the means eight readings each three plates. Only the displacements 
the alternate strong lines arising from the 4047 mercury line were mea- 
sured. The wave number shift the vibrational Raman line was obtained 


from excitation both the 4047 and the 4358 mercury lines. 


the time that these exposures were taken, arrangement had been made 
place intensity calibration marks the plates. However, after the measure- 
ments the positions the lines had been made was felt that the weak lines 
the rotational spectrum were sufficiently resolved and intense justify 
attempt measuring their intensity relative the strong lines. Therefore 
step weakener was made and calibrated electronically 4047 with the aid 
monochromator. The step weakener was placed the slit the spectro- 
graph and illuminated uniformly means ground glass screen and lamp. 
103a-0 plate the same batch used for the Raman spectrum was then 
exposed the spectrograph for hr. The calibration obtained was used 
determine the relative intensities the rotational lines the two plates. 
The densities the spectra were measured Moll microphotometer 
adapted with photo-multiplier cell and Brown recorder. 


Results and Discussion 
One the rotational Raman spectra obtained the way described 
above reproduced Fig. 2(a). 
The following method was used assign the correct values the rota- 
tional lines. The approximate formula given (see Herzberg (9)) for the rota- 


4 
— 
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The separation between successive lines 4B. the present case the position 
only every other line was measured that the separation between the suc- 
cessive lines given Column Table 8B. The observed value for this 
separation approximately 7.0 the other hand the separation 


TABLE 
DISPLACEMENTS AND INTENSITIES ROTATIONAL RAMAN LINES 


(vac) (vac) Relative Intensity 
observed calculated intensity ratio 
86. 86.11 2.5 
72.13 72.15 2.4 
58.17 58. 2.6 
44.10 44.09 3.3 
28. 30.00 2.1 
—22. —22.95 
—44.00 —44.09 3.4 
—51.07 —51.12 4.4 
—79. —79.14 2.8 
Mean 2.8 


Stokes line from the corresponding anti-Stokes line 8B(J 3/2). 
Therefore with the approximate value just determined the value any 
pair corresponding lines can determined. For example, from the two 
lines Av= 44.10 and 44.00 one obtains 3/2 88.10/7.0 


4078 


(a) 


891.85 


(b) 


Fic. print the Raman spectrum fluorine. 
(a) The rotational spectrum arising from the 4047 line. 
(b) The vibrational line 4535.7 arising from the 4358 line. 


: 
4358 
| 
a 
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12.6, that is, 11. this way the numbering all the other lines 
given (see Table I). 


the correction due centrifugal stretching included the formula for the 
displacements the rotational Raman lines becomes (see, e.g., Herzberg (9)) 


where approximately (2) 


Dividing Equation (1) 3/2 gives 


Therefore, plotting against 3/2)? should give straight line 


whose slope The experimental values are plotted against 
3/2 


the Stokes and the anti-Stokes lines the same value (Table I). 


3.550 
3.530 


3.510 


(J + %) 


° 200 400 600 800 


anti-Stokes lines given Table The slope the line was made equal 

seen that the accuracy the measurements not sufficient determine 
the slope the straight line with any accuracy. appeared preferable 
employ the theoretical value from Equation (2) using preliminary values 
for and and determine that straight line slope 0.000027 
which fits best the experimental points. This line shown Fig. The 
equation this line 


The values calculated this equation are given Table Column 


3.5310 0.000027 3/2)? (5) 


) 2 
° 
fe) 
2 
3.490 
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According Equation (3) the intercept this line with the ordinate axis gives 
Since 0.000020 this yields 


Using the fundamental constants adopted Herzberg (9) and 9.50227 


Hg 4046.6 A 


width 


Anti- Stokes Stokes 


50 


| 


Fic. trace the rotational Raman lines fluorine. The numbers 
the maxima are values. 


Fig. shows microphotometer trace the rotational Raman spectrum 
fluorine. The relative values the intensities the lines are given Table 
Column and the ratio the intensity line odd the mean the 
intensities two neighboring lines even given Column Owing 
the rough nature the intensity measurements the observed intensity ratio 
varies from 2.1 4.4 with mean 2.8. This value satisfactory agree- 
ment with the theoretical value 3.0 for nuclei spin Gale and Monk 

(6) well Aars (1) obtained values for the intensity ratio which varied from 
5.5:1. Since they were observing the emission spectrum fluorine, 
perturbations between the different excited electronic states may have been 
present. Nevertheless since the theoretical maximum value this ratio, 


they accepted the value for the nuclear spin Recently Gilliam, Ed- 


us 
wards, and Gordy (8) observed means microwave spectroscopy that the 
nuclear spin The evidence from the Raman spectrum 


agreement with this result. 


However, until now the statistics the nuclei had not been determined. 
From the fact that the odd numbered Raman lines are the strong ones and 
the assumption (which can hardly doubted) that the ground state the 
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molecule state must concluded that the nuclei follow Fermi 
statistics. This experimental result agreement with the general rule that 
nuclei odd mass number follow Fermi statistics. 


spectrogram showing the vibrational Raman line reproduced Fig. 
the present experimental arrangement the linear dispersion the plate was 
times higher than the arrangement used when the line was first observed 
(2). This resulted much higher accuracy. The value the displacement 
the vibrational line now found 


and compares well with the value 
previously reported (2). The first vibrational quantum its ground state 


Since Murphy and Vance (12) applied Badger’s rule (3) estimate the 
vibrational frequency, interesting ascertain how well the rule holds 
when the new values for the internuclear distance and vibrational frequency 
are used. The rule expressed 


where equilibrium force constant dynes per centimeter, 
0.68 
Since the values for these two constants equilibrium are not available 
present, the values used are those obtained the present experiment; that 
4.45 10° dynes/cm. 


Substituting and solving for gives 0.670 which agrees quite well with 
0.68 the value for required the rule. 
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EFFECT VARIABLE MASS THE ELECTRON 
THE SPACE-CHARGE LIMITED CURRENT DIODE! 


VISVANATHAN 


Abstract 


The change the current-potential distribution due the relativistic vari- 
ation the mass the electron has been calculated suitable series expansions 
the case plane parallel diode and has been shown considerable the 


case large power tubes. 
Introduction 


The space-charge limited current between the electrodes varies 
(neglecting initial velocities), where represents the potential difference and 
the distance between the anode and the cathode. Langmuir (2) deduced 
means simple dimension theory that dependence independent the 
geometry the system. This conclusion has been widely quoted the Three 
Halves Power Law. Experimentally, Chaffee has studied the current-potential 
characteristics various three-element vacuum tubes, connecting the grid and 
plate together order obtain diode system complicated geometry. 
has shown that under space-charge limitation the current can always repre- 
sented where near 3/2, over part the characteristic, though 
some cases the value must changed order fit the experimental 
data. The deviation from the Three Halves Power Law can due various 
causes, the most important one being the effect the variable mass the 


electron high potentials. 


The current-potential characteristic predicted Child’s law (1) does not 
take account the change the mass the electron due its increase 
velocity traverses the interelectrode space. This change, though very 
small ordinary vacuum tubes, may considerable large power tubes, 
where the potentials are high. shall calculate the magnitude the effect 
the case plane parallel diodes this paper. 


Plane Parallel Diode 


The motion the electrons governed the equation 


Manuscript received September 29, 1950. 
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This equation integrated subject the boundary conditions 


and 
(5) 


The integral can also written 


where 


The above integral can evaluated terms sum algebraic functions 
and elliptic functions the first, second, and third kinds. order get 
idea the effect and also for computational purposes integrate suitable 


series expansion. 


Assuming (usually true), expanding the integrand the Binomial 
theorem, and integrating term term, obtain 


The above series converges sufficiently rapidly that can computed 
for all values 0.8 within three significant figures. 


where the constant can evaluated from the first series Equation (7). 


m 
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Again expanding the Binomial theorem and integrating term term, 


From this equation, the value for can determined. For values 


which quite good for computation 1.1. When 1.1 the inte- 
grand can expanded decreasing powers and 


Thus see for large values clear that the electrons would move with 
uniform velocity. However, the range governed Equation (7) 
accompanied prohibitively high potentials for vacuum tube work. 
Now find the potential distribution throughout the space charge region, 
apply the energy law: 


Writing (4) 


(12) 


and seen above, only small values are applicable vacuum tubes. 
Hence keeping only the first two terms (7) may write 


como ot 3 13 


0c? 


then expanding the Binomial theorem and solving for keeping only the 
first two terms, 


V 
(14) 


Now substituting (14) (13) and solving for the current 


¢ 
= 
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This equation the same Child’s law except for the term the bracket 
which corrects for the mass variation. 


specific case for 100 kv., the deviation calculated with the 
6§e? V2 


term also obvious that one can neglect the initial velocities 


emission. similar calculation the case cylindrical diode contem- 
plated for future date. 
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THE STATIC DIELECTRIC CONSTANT DIPOLAR 


Abstract 


application general formula for the static dielectric constant 
made material having cubic arrangement dipolar molecules, each 
which has two equilibrium positions 180° apart and ordering forces the 
short range type which tend make nearest neighbors antiparallel. shown 
that such model cannot lead sharp transition dielectric constant unless 
changes lattice dimensions occur. Qualitative comparison with certain ex- 
perimental results made. 


Introduction 


Various discussions the static dielectric constant crystalline solids 
which the molecules can rotate have been given the literature. Thus 
(5) has calculated the dielectric constant for ketone the assumption that 
each dipole has two equilibrium positions and that the probability tran- 
sitions between these depends upon the applied field and upon dipolar inter- 
action between neighbors. Frank (4) has used similar two-position model but 
has assumed that the field acting upon given dipole proportional the 
mean resolved moment within the crystal. thus obtains second order 
transition using the Bragg—Williams approximation. The purpose the pre- 
sent paper give results calculation static dielectric constant using 
model the same type, but assuming (i) that dipolar molecules are arranged 
cubic structure, and (ii) that the more important forces between neigh- 
boring molecules are the short range type (nondipolar). The energy 
particular molecule result the action its neighbors taken into 
account, using the first approximation Bethe, that the effect local order 
upon the dielectric constant considered. 


Applicability the Model 


The importance short range forces dielectric materials has been dis- 
cussed Smyth (9) who has pointed out that some crystalline solids, glasses, 
etc., show dielectric behavior which evidently due dipole orientation, the 
restriction which dependent upon molecular shape, packing the lattice, 
and internal energy,—or, put another way, the restriction dipolar motion 
due largely short range forces between molecules. This restriction motion 
will normally such that each molecule will have several possible equilibrium 
positions which may, may not, have the same energy. The number such 
positions will depend the crystal structure and molecular shape. However, 
for some materials rather special crystal structure and molecular shape there 
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may exist two positions—having dipole directions differing 180°—-whose 
energies the completely disordered state are approximately equal, and con- 
siderably lower than the energy corresponding any other position the 
molecule. model based this type structure will used this paper. 
will also assumed that changes permanent moment the dipoles, 
electronic polarization, and lattice dimensions are very small over the tem- 
perature range considered. The model admittedly highly idealized one and 
there will probably materials conforming closely all these require- 
ments. However, for one two special substances may possible obtain 
reasonable agreement with experiment and the results should qualitatively 
applicable materials different crystal structure. 


Formula for the Static Dielectric Constant 


The method used calculating the static dielectric constant that given 
(6, and generalization Kirkwood’s method for dipolar 
liquids. Fréhlich’s method spherical region inside the dielectric medium 
considered, within which the dipoles are treated according the rules 
classical statistical mechanics. This spherical volume large compared 
region just big enough have the same dielectric properties macroscopic 
specimen. The region outside the sphere treated macroscopically con- 
tinuous medium. 


The spherical region then divided into unit cells, each which makes the 
same contribution the mean square moment the sphere. The static di- 
electric constant given 


where the optical refractive index, the electric moment the unit 
cell for given set (x) displacements its elementary charges, m*(x) the 
average moment spherical region volume inside the macroscopic 
specimen the set displacements unit cell the center (leading 
supposed persist for long time. the average value 
over all possible displacements (x) the unit cell and given 


where the probability finding the displacements (x), average 


over all possible states the surroundings outside the unit céll. No= N/V 
the number unit cells per unit volume. 


Formula (1) based the assumption that the electronic contributions 
may separated from those due rigid dipoles, etc., and treated 
macroscopic basis. This requires that the high-frequency dielectric constant 
not altered when the molecule turned from one position another, and 
further, that over the temperature range considered there change 
lattice spacing the material. The latter condition may rather stringent 
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one, since transition occurs which, practical cases, almost certain 
produce change crystal dimensions and may affect the electronic charge 
distribution considerably. This will discussed further when the application 
the model actual materials considered. 


Detailed Calculations the Dielectric Constant 


The type material considered this calculation one which dipoles that 
are nearest neighbors have opposite directions their lowest energy state, 
that the absolute zero temperature when the lattice completely ordered, 
the part the dielectric constant due rigid dipoles zero. the temper- 
ature increased, some the dipoles turn from their position lowest energy, 
and the numbers each type changing orientation must equal. thus 
possible consider each set dipoles forming one constituent mixture 
alloy which the quantities the two constituents are equal, and the 
problem calculating the change order with temperature the same that 
solved Bethe (1). Use will therefore made his results. 


case the lattice positions are divided into two types, and 
and dipole right (R) position points in, say, the 
direction, and wrong points the direction. The opposite 
true for dipoles the positions. The interaction energy pair 
neighbors neighbors zero, and pair neighbors has interaction 
energy 

The ratio the probabilities corresponding the two configurations the 
Boltzman factor exp The factor first approximation 
including long-distance order now called since the symbol has already 
been used for the dielectric constant. 


The sphere now divided into unit cells which contain one dipole each 
type, and the quantities Equation (1) calculated. picture one cell and 
its immediate surroundings for the case which each molecule has six nearest 
neighbors (simple cubic NaCl type) shown Fig. The cell dipoles are 
and When right its moment and when wrong 
The reverse true for 


Bun 


at 

@ 


Fic. cell (enclosed dotted line) and immediate surroundings simple cubic type) 
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The moments the unit cell for the four possible configurations are shown 
Table 


TABLE 
configuration moment configuration 


The probability corresponding the second configuration 
the probability that given (central) molecule and certain one its 
(first shell) neighbors This given 


n=0 


where the number neighbors the first shell (in addition 
for certain configuration that shell, 
the total number nearest neighbors 
the total probability that molecule and equal 
C(1 
factor (see below). 


The probability corresponding the third configuration similarly 


where the total probability that molecule and equal 
C(x 


Pwr == 


Since 


and since, from reference (1), can shown that Pwr 


should expected, since, the treatment consistent, the same result 
should obtained whether considered the ‘‘central’’ molecule. 
The remaining two probabilities may similarly determined and furnish 


There remains the problem evaluating for the configurations, 
and here approximation the moment the part the sphere outside 


1 = 
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the unit cell must obtained. doing this, care must taken treat equal 
numbers and molecules the same method. The following approx- 
imations will considered: 


Short range interaction between the unit cell and the remainder the 
spherical region neglected. Hence the moments all molecules outside the 
unit cell are considered straight probability basis. The average moment 
the molecules outside the unit cell is, therefore, 1)(r since 
the number unit cells the spherical region. The average moment 
the molecules similarly-(N 1)(r that the total moment 
the spherical region equal that the unit cell. 

Short range interaction between the unit cell and its immediate sur- 
roundings taken into account that the first layer molecules outside the 
unit cell Fig. included. The influence the outer molecules molecules 
factor for each wrong molecule this group,—that considered the 
central atom Bethe’s method. The molecules 10, therefore, belong 
‘‘second shell’’ around The factor cannot introduced for 
these molecules and the effect the outer layers them cannot considered 
this method. The influence this fact the accuracy the final result 
will considered later. 


method which would undoubtedly improvement that this 
section would one making use second approximation. The factor 
would then introduced for molecules 10, and the remainder 
would form the central group. This method considerably more complicated 
than the two described above, however, and will not attempted this paper. 
hoped that results its application may presented later. 


Calculation mm* for the second approximation 


The group first considered. The data for this group are shown 


Table 
TABLE 


Probability Moment 


The average moment is, therefore, 


(3) 


| 
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There are four such pairs, and, since they are independent each other and 
interact only with AB, their average moments may added. For molecules 
and the average moments are 


Hence 


The average moments are calculated exactly the same manner case 


(a) that 


Since the remaining two configurations and give zero moment for the 


unit cell, the over-all average given by: 


since 


Another approximation that Bragg and Williams (2, which long 
range order only considered. This treatment differs from those above that 
the probability given configuration the unit cell must now considered 
independent the interaction between its elements. The probability 
rw. Since the molecules outside the unit cell are also considered straight 
probability basis, their total moment zero. 


since and are now defined where the long range order 
Bragg and Williams. This quantity has been plotted accurately enough for 
the purpose this paper Nix and Shockley (8). becomes zero tran- 
sition temperature 1.2167 where the Bethe critical temperature 
which long range order disappears. 


The results for the three approximations may summarized follows: 
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where has the following values for the simple cubic case: 


for the first approximation taking account short range order. 


* 


for the second approximation taking account short range order. (Note that 
this still refers Bethe’s first approximation but takes more accurate ac- 
count the interaction between unit cell and the remainder the spherical 


region. 


Limits Error. the Approximation 


Before discussing the results obtained with these approximations im- 
portant consider the directions and magnitudes the errors involved. The 
probabilities used for the various configurations the unit cell will con- 
sidered first. For the first two cases these are calculated directly from the 
results Bethe’s first approximation and suffer from the same limitations. 
These have been discussed elsewhere (8) and sufficient point out here 
that the approximation generally considered reasonably accurate 
one. For the Bragg—Williams case the complete omission short range order 
makes the error more serious and has the effect raising the critical temper- 
ature 22%. 


The errors involved the calculation the moment the whole sphere 
when and are fixed now considered. Since interaction between molecules 
assumed make them unlike, the unit cell, for the case will 
the first approximation this effect neglected, that when 
the molecules are less (more positive) and the molecules 
less (also more positive) than they should be. Hence, since the moment 
the unit cell positive this case, the value too 
large. Similar reasoning shows that when the value 
used for more negative than should that again too 
large. should, therefore, expect the first approximation give value 

that too large, and, since all interaction between the unit cell 
and its surroundings has been neglected, this error may considerable. This 
same result refers also the approximation, since the method 
calculating the same. 


For the second approximation the molecules have been considered 
with accuracy corresponding that Bethe’s first approximation, that 
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the influence the outer molecules them have been taken into account. 
For the molecules 10, however, the influence the outer shells has 
been neglected. This influence will the same direction for all configura- 
tions the unit cell, and were equally important for 
for the net result mm* would zero. should, how- 
ever, expect that this influence will more important when opposing 
the molecule than when aiding, that there slight net effect tending 


for W). This net effect has been neglected the 


approximation, that the value smaller than should be, although 
probably quite close the true value. Thus, while the first approximation 
and the methods both give results that are considerably above the 
true value, the second approximation gives one that slightly below it. 


Results for the Simple Cubic (NaCl Type, 


The curves vs. are plotted for the three cases Fig. 
and the second approximation using the Bethe factor replotted larger 
scale Fig. These plots all show discontinuity slope (at 


2.0 


z 
0 1.0 2.0 3.0 
v/ To 
(a) Bethe first approximation. first outside layer with unit cell not 
considered. 
(b) Bethe first approximation. first outside layer with unit cell 
included. 


(c) approximation. 


but for the most accurate approximation (Curve this very 
slight. Furthermore, the change slope the transition temperature for 
the opposite direction those for and and this may probably taken 
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confirmation the fact that the approximations are opposite sides 
the true value. The position the latter perhaps very roughly similar 
that shown dashed. also seems probable that the true curve 
e € 
course, due the fact that long range effects apply equally the two types 
dipole and, since for the model chosen these have opposite directions, they 


n 


cancel each other and the net effect the dielectric constant zero. Thus, 
for this type model (subject the restrictions described the second 
section above) there should sharp transition the curve dielectric 
constant. The presence sharp transition, result approximate 
calculation, merely indicates that part the short range interaction has not 
been taken into account. 


Comparison with Experimental Results 


properties may probably interpreted qualitatively terms the model 
this paper. Accurate comparison hardly justified since, view the shape 
the molecule, one may question whether the crystal structure will 
simple cubic the type here discussed and whether there will only two 
equilibrium positions 180° apart. However reasonable expect that the 
general forms the calculated and experimental curves will similar. 


The curve obtained White, Biggs, and Morgan (11) for the static di- 
electric constant this material over range temperatures shown 
Fig. These authors infer from the shape the molecule that rotates about 
one axis and, using model proposed White (10), assume that there are two 


5.0 
1.0 2.0 3.0 4.0 
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potential minima which differ value and are 72° apart. They show that the 
curve polarization calculated using this model fits reasonably well. 


determined from the curve Fig. proportional the quantity plotted 
Fig. and proper choice approx.) the experimental data 
can made produce curve which has the general form Curve Fig. 
Detailed comparison would probably pointless however since, addition 


Using reasonable value the quantity can readily 


STATIC DIELECTRIC CONSTANT € 
o 


-80 160 240 
TEMPERATURE °C. 


Fic. dielectric constant (60 sec. charging time) 6-trichloro- 
benzene determined White, Biggs, and Morgan (11). 


the factors mentioned above, neither the exact form the true curve 
9 
(shown roughly the dashed curve Fig. 2), nor 
€ 
known. possible that the latter could determined specific heat 
measurements. Such measurements would have quite accurate however, 
the change configurational specific heat the transition temperature 
becomes small fraction the total specific heat for fairly large molecules 
the type here considered. Further difficulties might arise the rate approach 
the equilibrium state were very low. 


The importance this rate approach equilibrium should especially 
emphasized. The order-disorder theory based statistical mechanical con- 
cepts which apply only the system equilibrium. many materials, 
even temperatures near the establishment equilibrium may require 
quite long periods time. This clearly brought out measurements 
these apparent that, for temperatures below 170°K. (which, this model 
applies, approximately even ballistic measurements using charging 
time sec. may error. Relaxation times the same order presumably 
apply thermal measurements. make the usual assumption that the 
relaxation time exponential function temperature the material may 
become essentially temperatures not very far below the transition 
temperature, and may not reach very highly ordered state. This effect may 
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reduce the measured value configurational specific heat considerably. 
any event much time possible should allowed for the material reach 
equilibrium low temperatures. 


Conclusions 


From the comparison the above theoretical results with experimental 
data may draw some general conclusions concerning dielectric transitions 
crystals whose structure such that neighboring molecules tend orient 
themselves with their dipoles antiparallel directions. material this 
type shows sharp transition, must connected with expansion the 
lattice which allows sudden increase freedom the dipoles (or the whole 
molecule) with, usually, appreciable changes electronic polarization. This 
expansion may accompanied change completely different type 
crystal structure with different number nearest neighbors given 
molecule. accurate theoretical calculation the variation dielectric 
constant transition this type can made only account can taken 
most these factors, and practically every analysis therefore becomes 
special case dealing with specific material. 
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ROTATIONAL ANALYSIS THE BAND SYSTEM 
THE MOLECULE! 


ZEEMAN? 


Abstract 


The band heads the spectrum discovered 1932 Fowler and Bakker 
have been remeasured plates taken large quartz spectrograph and first 
order grating plates. The (0,0) and (0,1) bands have been photographed 
the 6th, 8th, and 9th orders concave grating, and rotational analysis 
these bands has been carried out. For the ground state was found that 


the value was obtained. 


Introduction 


1932 Fowler and Bakker (2) discovered new spectrum discharge 
through sulphur vapor containing nitrogen impurity. consisted two 
band systems bearing striking resemblance the and systems NO, 
both regards structure and wave length. view experimental and theo- 
retical considerations they attributed these bands the molecule. Thus 
far rotational analysis these bands has been made, and the object 
the present paper submit such analysis the bands. 


Experimental 


Since the and band spectra lie nearly the same spectral region 
essential eliminate the bands completely from the light source 
order photograph the band systems successfully. Fowler and Bakker 
used commercial nitrogen, dried passage over phosphorus pentoxide, and 
produced the bands passing gentle stream this dried nitrogen 
through ordinary end-on discharge tube while the sulphur reservoir was 
gently heated with Bunsen flame. the present research pure nitrogen pre- 
pared from pure NaN; and this gas was circulated mm. 
mercury pressure through ordinary end-on quartz discharge tube means 
two-stage mercury diffusion pump. Sulphur was contained vertical 
reservoir attached side tube protruding horizontally from the capillary 
the discharge tube. was found essential use vacuum distilled sulphur 
eliminate bands which appeared preliminary plates where 
ordinary flowers sulphur were used. The partial pressure the the 
discharge tube was adjusted until the second positive nitrogen bands and the 
visible bands were about equal intensity. This condition could kept 
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constant for hours means cylindrical electrical heater which surrounded 
the side tube. 


Molybdenum electrodes were used instead aluminum because the latter 
metal combined with nitrogen higher temperatures causing the pressure 
the tube decrease unduly. The capillary the discharge tube was cm. 
long, and had mm. with wall thickness 2mm. k.v.a. trans- 
former giving 6600 the secondary circuit was used pass the discharge 
through the tube. The primary current was kept constant amp., using 
220 supply. 


was found impossible eliminate the persistent bands completely 
from discharge through alone even though this gas was circulated over 
heated copper turnings. However, when sulphur vapor was added the 
bands did disappear entirely. The final exposures were made type 
Eastman spectroscopic plates using the sixth, eighth, and ninth orders 
grating Paschen mounting. The exposure time was about hr. 


order photograph these high orders was necessary eliminate all 
wave lengths above about 2500A from the light entering the slit. This could 
not done means material filter, and the cylindrical lens method de- 
scribed Herzberg and Curry (5) was employed. However, eliminate 
completely the very strong and band systems was found necessary 
introduce adjustable preliminary slit near the window the discharge tube. 


Vibrational Analysis 


The bands originate from transition quite similar 
the bands NO; they are degraded the violet. The bands result 
from transition, and are degraded the red. 


The heads these bands have been remeasured first order grating 
plates and plates taken large Hilger spectrograph. Essentially the same 
Deslandres diagrams given Fowler and Bakker were obtained with the 
new measurements. The heads the bands can represented the fol- 
lowing 


The heads the bands can represented follows:— 


Vaead = 39882.4 (953.350’ =< 8.650'?) 


' 
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Rotational Analysis 


Fig. spectrograms the (0, and (0, bands are reproduced. 
some distance from the heads the branches were resolved into different series 
lines which could followed towards the heads using constant second 
differences. studying the relative intensities these series lines, the 
nature the branches could established. the sub-bands the 
and branches, and the sub-bands the and branches 
were found the strongest both the bands analyzed. This the same 
result was obtained Guillery (3) the analysis the bands. The 
(0, band much stronger than the (0, band, and higher degree 
accuracy was obtained the analysis the former band. The plates were 
measured Abbé comparator, and the wave numbers computed using 
the method described Almy and Horsfall (1). 


The ground state intermediate between case (a) and case 
with fairly large value for the coupling constant For such cases Hill and 
Van Vleck (6) bave obtained the following formula for the spectral terms 
involved 


For large [case expansion Equation (5) gives:— 


for regular doublet splitting with the sign for the component and the 
sign for the component state. 


The formula for the rotational terms can further simplified omit 
the terms independent and also introduce the value 
denoted (7), that for doublet states 


substitution Equation (6) and dropping the subscript the following 
formulas are obtained for the rotational terms the two sub-levels the 


the present investigation the bands splitting the upper state 
was observed, and thus the expression 


of 


(0,0) 


grating ina Paschen mounting. 
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was sufficient for the rotational terms this state. Hence omitting the term 
and the subscripts have for the state, 


3% =F; 
+F, 
2 +F, 
(4) 
(3) 
| 
(2) c+ y, 
(2) 
2 


Fic. energy level diagram for the first rotational lines band, 
showing the transitions for the determination 


The broken lines represent the satellite branches for case (b). The branches and are 
forbidden case 

The spin-doublet splitting the lower state much larger than shown, and the spin-doublet 
splitting the upper state was not observed the present case. The A-doubling splitting the 
state negligible comparison with that the state. 


The designation the extreme right-hand branch should instead 


| 
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According Equations (8), (9), and (10) the expression for the term values 
involve for states, and for states. However, the levels state 
can formally numbered even for small values. this paper the 
combination relations dealing with the state are given terms and 
those the state terms Hence facilitate the evaluation the 
combination differences and the rotational constants, the band lines were 

Tables and the wave numbers the band lines the (0, and (0, 
bands are given. Table the combination differences the upper states 
the (0, and (0, bands are given. These quantities were determined 
from the following 

(12) 


The combination differences for the lower states were obtained using the 
following relations:— 


and 


Because the A-type doubling that occurs for all states with the lines 
the and branches under these circumstances always end somewhat 
different states from the lines the branches. Hence the so-called ‘‘com- 
bination occurs the and values of, for example, 
transition, that these quantities cannot used determine the 
rotational constants accurately. However, the values were formed 
shown below, order obtain indication the magnitude the A-type 
doubling. knowledge thereof enables decide whether the splitting the 
state normal inverted. According theory (4, 263) the A-type 
doubling for the sub-band considerably smaller than that for the 
sub-band. the present investigation, the combination defect was found 
the order the experimental error for the sub-bands with smaller wave 
numbers while was considerably greater for the other sub-bands. This indi- 
cated that the spin-doublet splitting the *II state normal for NS. The fol- 
lowing relations for example were used evaluate the combination defect 


for the sub-band 


for the sub-band. 


—_ 
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TABLE 


Kk" 


42.42 
123 41.41 
143 
15} 39.12 
163 38.68 
38.18 
38.03 
37.86 
38.18 
243 39.12 
40.50 
41.99 
43.10 
314 45.61 
323 46.94 
334 
50.2 
52.04 
57.72 
48} 


59.11 67.77 
59.11 74.40 
59.51 76:13 
82.71 
85.14 
70.7 
76.47 16.15 
78.63 
83.11 
63.94 
68.96 
20.51 84.99 
90.47 
28.51 
41.42 13.69 
45.64 19.81 
51.91 


41974. 
76. 
80. 
84. 
87. 
90. 
94. 
97. 

06. 


74.39 81.37 
72.42 81.03 
70.98 80.85 
80.56 
67.77 80.85 
81.37 
81.92 
64.68 82.78 
83.54 
62.54 85.98 
61.96 
61.64 
61.64 
61.64 90.99 
61.64 93.25 
61.80 94.59 
61.96 
62.54 98.79 
63.75 
64.68 05.73 
65.88 
66.98 11.09 
13.91 
69.39 16.87 
70.98 20.07 
72.7 
74.39 
76.26 30.22 
80.56 37.71 
82.19 41.47 
85.20 45.55 
50.10 
89.76 54.22 
92.79 58.64 
63.11 
72.47 
77.24 
82.65 
98.41 
15.64 
21.45 
46.21 


57.60 
61.96 
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99.05 
11.37 
17.56 
24.07 
30.41 
37.25 
44.32 
51.06 
58.00 
65.77 
72.42 
80.67 


87.81 
89.32 
90.80 
O6 
12! 
74.96 
283 
71.48 
313 
92.08 
97.54 
373 08.78 
14.53 
403 26.58 
41} 32.90 
42} 
51.97 
58.92 
65.77 
72.42 
79.47 
503] 
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TABLE 


| | 


| | 


| 
K” | J" | Rs | Py R, | 

103 46.77 67.98 86.03) 05.89 27.10 
123 44.10 64.14 85.93 133 66.11 35.46 
143 42.76 66.01 65.00 89.31 43.21 
42.32 66.80 163 90.43 18.38 48.19 


to 
for) 
or 
ou 


TABLE 


bas | 
_ 


523 72.09 

81.50 
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These relations, well those already mentioned for the calculation 
and can easily checked the energy level diagram 
Fig. 

evaluate for the various values for the state, use was made 
the following method. substituting the expressions (8) (9) for 


straight line obtained. From the intercepts these lines with 
the ordinate axis the values and can obtained, and from the slope 
the corresponding values. Finally obtained taking the mean 
value and BY? suggested Equation (7). 


The constants and can obtained from the combination differences 
manner similar that described for the state. 


this paper semifinal values and were determined this method, 
and the final values and The determined value, 1.22 


agrees well with the expected value 1.23 calculated 
3 
from the formula using the new and values. However, the 


2 


accuracy the data for the (0, band not high enough determine 
reliable value. 


Using the value thus obtained, and approximate value de- 
against values was constructed obtain accurate value 
order obtain the best possible values and the differences 


However, the ordinary relations (4, 186) cannot used for 
transition, because additional terms now appear these formulas. Fortun- 
ately the differences B,’) can still obtained appropriate choice 
branches indicated the following theory. 


From Equations (8), (9), and (11) follows that for transition 
the present type, 


(19) 
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Thus 


and were determined for the (0, and (0, bands. However, 


develop the following additional formulas. 
From the above 
Similarly can proved that 


(22) 


Hence plotting the average the obtained Equation (22) 
determined for the bands analyzed. 


means the procedure explained above, the following values were finally 
obtained: 
(0,0) (0,1) 


and 


(0,0) (0,1). 
These values can checked the following way. For the same levels the 
(upper) state, follows from the Equations (19) that 
(0,0) (0,1) (0,0) 
and 
(0,0) (0,0) 
(0,0) (0, 
can obtained. The following values were determined this 


values agree well with the method already described. 
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Evaluation the Coupling Constant 
From Equation (6) follows that for state) 


2B? 


The differences Equation (25) are identical with the quan- 
tities Ri(J)— Q2(J), already mentioned above (Equation (22)). Hence 
comparing Equations (25) and (22) 
3B; 
(1) (2) 
(26) 
this way the coupling constant was determined. example this 
graphical determination the graph for the (0, band given Fig. 


Table summary the rotational constants the bands given. 


TABLE 
SUMMARY THE ROTATIONAL CONSTANTS THE BANDS 


0.75 200.75 400.75 600.75 800.75 1000.75 1200.75 1400.75 1600.75 1800,75 2000.75 2200.75 2400.75 2600.75 
J 


evaluate the coupling constant explained the text. 


237 
236 
235-- 
233 
232 
230 
227 
226 
225 
224 . 
223 
222 
fi 
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For the ground state, finally have the following constants* for the 


molecule; 


Also, the folowing values were obtained for the constants listed 


For the (0, band: 43387.35, and 43165.88 
For the (0, band: 42183.24, and 41961.71 


The coupling constant 223.03 
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Abstract 


search has been made for the reaction (n, have been unable 
find any activity ascribable the energy greater than Mev. 
and the half life lies between sec. and min., the capture cross section for the 
above reaction less than microbarn. 


Introduction 


The carbon isotopes found nature have mass numbers and 13. Carbon 
and carbon are short-lived positon emitters with half lives sec. 
and 20.5 min. respectively. the neutron-excess side the stability curve 
lies long-lived negaton emitter with half life variously estimated 
between 5100 and 7200 years. Until recently mention has been made the 
literature carbon mass 15. 


From the mass calculated Bethe (1) for can 
seen that 10.8 Mev. difference exists between and the stable end product 
mass 15.0049. Using the data assembled Konopinski (5), one would 
expect roughly the following all the energy were expended transition: 
(a) were allowed transition, the half life would about 0.04 sec., 
(b) first order forbidden, about sec., (c) second order forbidden, about 
min. Using recent value for the mass Tollestrup, Fowler, and 
Lauritsen (8) 14.007733 amu and the mass for calculated Bethe, one 
would expect that the reaction would occur with the release about 
0.2 Mev. 


Recently, Hudspeth, Swann, and Heydenburg (3) have reported the pro- 
duction the reaction (d, They bombarded and 
active barium carbonate with deuterons from Van Graaff generator and 
obtained, difference, emitter with half life 2.4 sec. and energy 
8.8 Mev. which they ascribed the transition the ground 
would then exothermic about Mev., instead 0.2 Mev. 


Hudspeth a/. not discuss chemical purification the barium carbonate 
which was used their experiment. Unless chemical identity was established 
between the active and inactive barium carbonate, the activity may have been 
due impurity the active material, since the processing the two mater- 
ials was surely different. 


connection with general program determining cross sections un- 
stable nuclides and finding new nuclides, have made search for The 
Manuscript received original form August 1950, and, revised, December 12, 1950 


Contribution from the Chemistry Branch, Atomic Energy Project, National Research 
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nuclear reaction looked for was y)C and the region half lives inves- 
tigated thus far has been sec. min. 


have been unable find any activity which might ascribed 
These results have previously been reported very briefly Letter the 
Editor the Physical Review (9). 


Experimental 
(a) Preparation Material 


The active carbon sample was received solution sodium carbonate 
containing atom The carbon both active and inactive sodium car- 
bonate solutions was reduced elemental carbon using the same technique 
and reagents. This assured the chemical identity the two samples, which 
necessary prerequisite the isotopic assignment made without chem- 
ical purification after irradiation. 


Barium carbonate was precipitated from the sodium carbonate solution, 
centrifuged, washed with water, then with acetone and gently dried. The 
barium carbonate was then mixed with lead chloride (2) and the mixture 
heated vacuo 350°C. where the lead carbonate which formed the dry 
reaction decomposes with evolution carbon dioxide. The carbon dioxide was 
then frozen down and purified repeated distillation. 


The carbon dioxide was reduced elemental carbon (6) the following 
manner. The center in. diameter iron pipe about in. long was filled 
with 10-fold excess C.P. magnesium turnings. One end this pipe was 
connected the vacuum system. The other end was capped and the system 
was made air-tight using Khotinsky cement and black wax. The reaction 
tube was evacuated pressure mm. and the tube then filled 
pressure cm. with carbon dioxide. The center section the pipe was 
heated about 650°C., the melting point magnesium, which temperature 
magnesium reacts with carbon dioxide, reducing carbon. The remainder 
the carbon dioxide was converted carbon allowing successive small 
amounts carbon dioxide into the pipe. 


When the reaction was complete the solids were removed from the pipe and 
treated with water and hydrochloric acid, gradually increasing the hydro- 
chloric acid concentration until was The carbon was leached for about 
one hour this manner, then centrifuged. After second treatment with 
hydrochloric acid for hour was centrifuged, washed with water, then 
acetone, and finally dried vacuo. 


Procedure 


Irradiation the rapid action rabbit the NRX pile requires the encasing 
the sample capsulelike containers. Since the half-life expected was short, 
comparison with the time required remove the sample from the capsule 
the activity the irradiated material had measured the capsule. 
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‘Spec-pure’ iron was chosen suitable material for the capsule construction 
since the capture cross section the iron low and all resultant activities have 
reasonably long half lives. ensure chemical identity, all capsules were 
fabricated from the same rod. 


special flask, shown Fig. was used for receiving and immediately 
counting the sample. The flask was cylindrical and made lead provide 
adequate shielding. The sample the capsule drawn above the posi- 
tioning plunger suction applied the valve locating plunger 


A- Top of Pile 

Shielding 

C- Vacuum Valve 

D- Locating Plunger 
E- Positioning Plunger 
F- Proportional Counter 
G- Absorber Frame 
H-Pneumatic Tube 


Fic. Castle for receiving and counting samples. 


which normally held position spring. When this plunger depressed, 
its length travel will indicate the presence absence the capsule. pro- 
portional counter with thin glass window (30 mgm. per filled total 
pressure cm. 20°C. with argon (67 cm.) and carbon dioxide cm.) was 
placed position shown. When position the capsule directly opposite 
the thin window, with intervening shielding. absorber frame placed 
between counter and capsule that the same apparatus may used for 
obtaining absorption curves and determining the maximum energy the 
emitters found. 


Suction <—Bd 
Py 
- 
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The entire assembly was placed over the pneumatic tube which leads 
high flux position neutrons per cm.? per sec.) the pile. 


The procedure followed was have the capsule held above the positioning 
plunger suction alone. When valve was closed the capsule dropped the 
bottom tube After the desired irradiation, the capsule was again sucked 
into counting position above reopening 


(c) Counter Assembly 


The counter assembly consisted preamplifier, linear amplifier, high vol- 
tage (3000 v.) set, scale 1000, and Recorder. The speed 
motion the paper the recorder was determined using 60-cycle line 
frequency. Each division (0.75 in.) corresponded 3.7 sec. the speed used. 
This time interval could further divided into smaller parts using dividers. 


The linearity the proportional counter assembly was checked and 
losses greater than were found counting rate 200,000 per minute. 


(d) 


The various irradiations made are shown Table 


TABLE 
IRRADIATION DATA 


Time 
Weight, gm. irradiation, 
sec. 


Results Material 
shown irradiated 


0.0013 


Results 


The results the irradiations are shown Figs. 2-5 incl. Fig. the 
decay curve obtained from the irradiation block the ‘spec-pure’ iron 
used capsule material shown. least three activities are present with 
half lives 4.4 0.2 sec., 26.5 0.3 sec., and 12.9 0.1 min. The limits 
shown are those internal consistency. These activities are probably due 
and mixture and The first four are produced 
(n, and 2m) reactions impurities known present, the fifth 
(n, 2n) reaction The activities found irradiation the inactive 
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and active carbon samples are identical those found the iron itself, 
shown Figs. and additional activity due the introduction 


was found. 


3 

2 


4 8 16 20 24 2 32 3% 7 4 4 #5 120 
TIME (1 UNIT = 3.7 SECONDS) 


COUNTS PER UNIT TIME 


wo OY 


Fic. Decay curve irradiated iron cylinder. 


COUNTS PER UNIT TIME 


4 8 12 16 200 28 32 40 4 Il2 "6 120 
TIME (1 UNIT = 3.7 SECONDS) 


Fic. curve normal carbon irradiated iron capsule. 


Since took approximately sec. for the sample leave the pile and the 
irradiation times were sec. duration, activity which had half life 
the sec. min. range would have been detected. 


3 
2 
26 SEC. 
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set upper limit the cross-section for thermal neutron capture, 
Li? (as lithium carbonate) was irradiated. value millibarns (4) was 
taken the cross section for the reaction The amount 


COUNTS PER UNIT TIME 


> 


44 


6 (20 
TIME (1 UNIT= 3.7 SEC.) 


Fic. Decay curve carbon containing irradiated iron capsule. 


COUNTS PER UNIT TIME 


2 24 2B 6 
TIME (1 UNIT = 3.7 SECONDS) 


Fic. Decay curve lithium carbonate irradiated iron capsule. 


irradiated (1.2 gm.) corresponds equal absorption neutrons for 
the weight used, assuming cross section microbarns for the latter. 
The 0.88 sec. 12.7 Mev. activity was readily detected shown Fig. 


SEC. 
107 
8 
4.7 SEC. 
0.98 
SEC. 
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eliminate the possibility nonresponse the counting equipment 
lower energies, lithium fluoride was irradiated. Both the and the sec. 
Mev. were found the amounts expected. 


Even though the results are reproducible better than 10%, was assumed 
that the minimum increase activity due which could have been de- 
tected was 20%. Using this reference, the cross section for capture 
was then calculated for various half lives the range studied. was 
assumed that the energy was greater than Mev. that significant 
absorption did not occur. 


For half life sec. the upper limit for the capture cross section 
0.05 microbarns. This increases 0.6 microbarns for half life for min. 


has eight neutrons the nucleus and thus one the magic number 
nuclei suggested Mayer (7) and others. Although one might expect the 
cross section low for this reason, yet our upper limit less than micro- 
barn times smaller than the smallest cross section reported for any nuclide, 
namely for 
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SUR CORRELATION ENTRE LES DISCONTINUITES DANS 
VALLEE DES MASSES LES ENVELOPPES NUCLEAIRES. 


Par BRONIEWSKI 


Sommaire 


montre que les discontinuités qui produisent dans variation régulitre 
vallée des masses fonction rattachent approximativement, dans 
partie tableau périodique, phénoméne saturation des en- 
veloppes protoniques dans partie, saturation des enveloppes 
neutroniques dont succession été établi par Mayer. 

L’hypothése d’une discontinuité dans région des éléments 
transuraniens est discutée. Une telle discontinuité aurait pour conséquence 
modifier les prédictions actuelles faites sur les energies des éléments trans- 
curiens. 


Introduction 


est bien connu que représente toutes les variétés nucléaires naturelles 
dans plan rectangulaire Neutron-Proton, obtient une multitude 
points compris dans une bande assez étroite variant premiére approximation 
d’une continue avec nombre masse 


sait aussi que les noyaux légers qui correspondent des valeurs 
ont tendance grouper autour bissectrice des axes mais que, par 
suite répulsion Coulombienne, tous les isotopes plus lourds 
plus plus cette bissectrice, leur nombre neutrons excédant celui des 
protons. Ainsi pour nombre isotopique A-2Z N-Z) atteint 
valeur 54. 


II. Etude des discontinuités dans vallée des masses 


Nous allons examiner maintenant plus détail vallée des 
masses dans plan diagramme qui semble faciliter cette étude est 
représenté sur Fig. fait figurer seulement les isotopes par plus 
léger plus lourd tous les éléments pair les points figuratifs ont 
été reliés par les lignes épaisses dans but faire resortir les limites sta- 
bilité dans chaque région. Avec Guggenheimer (6), nous appelerons ces limites 
les limites extérieures D’autre part, les noyaux impair qui, 
n’étant représentés que par par deux isotopes stables, sont particuliére- 
ment interessant car ils marquent trés sensiblement fond vallée des 
masses (minimum des paraboles Bohr-Wheeler) ont aussi été représentés 
reliés entre eux par des lignes pointillé (2). tel diagramme rend compte 
clairement double tendance qui préside des noyaux stables. 
sait effet qu’approximativement dans premiére moitiée vallée des 
masses les noyaux stables moyenne par adjonction d’un neutron 
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pour chaque proton ajouté, tandis faut deux neutrons pour chaque proton 
additionnel vue une configuration comparablement stable 
point vue radioactivité dans derniére partie tableau périodique. 


Fic. Diagramme indiquant les limites stabilité vallée masse leur liaison avec 
les enveloppes nucléaires dans schéma Mayer. 


Lorsque premiére ces tendances intervient voit, sur Fig. que 
les noyaux comparables entre eux point vue radioactivité (c.a.d., ceux 
qui impairs sont situés fond vallée masse méme les noyaux 
qui marquent les limites stabilité) placent trés exactement sur 
les droites pente lorsque c’est deuxiéme tendance qui devient prépon- 
dérante, ils placent sur les droites pente C’est d’ailleurs considérant 
regularité remarquable dans des noyaux lourds trés lourds, 
les moins riches neutrons, (tombant sur les lignes inférieures 
Fig. que Duckworth collaborateurs (4) ont découvrir nouvel isotope 
stable platine, probable plusieurs autres isotopes rares 
qui sont suivis points d’interrogation sur Fig. 


L’examen Fig. montre cependant que régime pente 
régime pente d’une continue qu’a une portion, assez limité 
vallée des masses correspondant environ dix éléments consécutifs. 
preuve est fournie par les brusques déplacements des limites extérieures 
stabilité qui produisent plusieurs endroits tableau qui seront 
étudiés détail par suite. Notons dés maintenant que ces déplacements 


32 SS 58 126 88 
sé 28 452 
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produisent pas méme maniére pour les deux régimes stabilité con- 
siderés. voit effet qu’avec systéme coordonnées les droites 
pente déplacent plusieurs reprises vers haut (en suivant sur 
vallée des masses sens des croissants), c’est dire, que discontinuité 
apparait pour une certaine valeur nombre atomique contraire 
régime pente régit stabilité nucléaire, vallée des masses décale vers 
droite, c’est dire que discontinuité produit pour une certaine valeur 


L’interpretation plus simple cette situation pourrait étre que ‘‘le ré- 
gime pente appliqué sans discontinuité une portion étendue vallée 
des masses, conduirait des noyaux trop pauvres neutrons ayant par 
conséquent tendance subir radioactivité Les réadjustements font 
donc dans sens qui accroit plusieurs vallée 
nombre isotopique Par contre, semble que les noyaux qui par 
deux neutrons par proton arrivent contenir aprés plusieurs 
semblables plus neutrons qu’il n’est requis par les diverses con- 
ditions stabilité nucléaire. d’autres points critiques des 
réadjustements produisent dans but cette fois diminuer nombre 
isotopique. 


Passons maintenant détermination exacte ces points critiques 
variation reguliére vallée des masses trouve interrompue. 


Ellis (5) rattaché des discontinuités aux nombres 20, 
neutrons qui correspondent aux ‘‘nombres bien connus aussi aux 
nombres 32, 100 neutrons. Ellis n’a pas cru devoir associer les 
des protons aux discontinuités dans vallée des masses. nous 
semble contraire qu’il faille ces discontinuités fois aux nombres 
magiques neutrons aux nombres magiques protons ayant soin toute- 
fois bien préciser les régions tableau periodique considére. 


produisant saturation des enveloppes protoniques 


Considérons d’abord vallée des masses dans région des élé- 
ments trés légers. sait que proton sont les deux seuls noyaux 
stables ayant nombre isotopique négatif, égal aussi ils placent sur 
portion droite pente située dessous bissectrice des axes 
Ensuite, jusqu’a l’oxygéne, existe régime les noyaux situent sur 
deux droites pente droite inférieure correspondant bissectrice des 
limite supérieure stabilité déplace encore vers haut d’une unité, les 
noyaux situés sur cette droite ayant nombre isotopique égale d’autres 
termes, vallée des masses subit deux premiers déplacements vers haut pour 
les deux premiers nombres magiques protons, Ces nombres n’ont pas 
été mentionnés par Ellis, effectivement aucun déplacement similaire n’est 
observable pour ces deux nombres magiques neutrons. 


continuant suivre vallée des masses dans sens des croissants 
voit produit d’abord, avec protons, déplacement limite 
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supérieure stabilité vers haut, deux unités ensuite, pour 20, 
autre déplacement vers haut méme importance mais intéressant, cette 
fois, les deux limites extérieures stabilité. 

peut essayer rattacher premiére discontinuité nombre magique 
neutrons comme fait Ellis (il faudrait alors admettre que déplace- 
ment fait gauche sur horizontal 20*) comme nous pré- 
férons fait que nombre correspond saturation des deux premiéres 
sous-enveloppes protoniques (comptant états possibles) s’édifiant 
vallée des masses produit pour neutrons deuxiéme discontinuité 
signalée ci-dessus doit incontestablement rattachée nombre magique 
protons. 

les limites extérieures stabilité obéissent encore strictement 
pente discontinuité qui produit pour est double nature; 
part limite supérieure stabilité déplace vers haut deux 
tendance des noyaux stables suivant mode que nous avons 
convenu d’appeler ‘‘régime pente est interessant remarquer qu’avec 
protons sous-enveloppe qui s’édifie partir calcium 20) com- 
prenant états trouve saturée. existe méme des indications expérimen- 
tales prouvant que configuration nucléons correspond arrange- 
ment nucléons stable (3). Notons aussi n’y pas 
discontinuité dans vallée des masses pour neutrons. 


pente continue que pour quelques éléments partir 
viron 46. partir cette valeur pour tous les éléments suivants 
rencontre plus que tendance d’adjonction deux neutrons par proton. 
peut noter qu’avec sous-enveloppe comportant états qui 
aprés nickel 28) trouve saturée. Par contre changement pente 
qui produit pour peut étre rattaché aucun effet saturation 
nucléaire tout moins dans schéma Mayer (7). 


Fig. montre que limite supérieure stabilité avec 38; 
fait doit ¢tre relié des isotopes strontium contenant plus 
neutrons. Cette perturbation dans vallée des masses produisant pour 
neutrons qui été signalée par Ellis, apparait cependant comme effet 
purement local changement pente translation générale 
fond vallée comme témoigne, dans cette région, l’examen ligne 
pointillée reliant les éléments impairs. 

discontinuité suivante dans vallée des masses, derniére produire 
suivant axe vertical, c’est dire, pour nombre donnée protons, cor- 


vers haut limite supérieure stabilité apparait déja pour 16. 

Nous sommes probablement ici dans une région tableau les limites stabilité trouvent 
influencées fois par une configuration stable protons neutrons. 


‘ 
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respond nombre magique protons. Fig. montre que cette dis- 
continuité consiste une extension vers haut deux unités limite 
supérieure stabilité. considération ligne stabilité maximum 
indiquée par les noyaux impairs montre que cet effet est aussi accompagné 
d’une translation analogue fond vallée des masses. Nous venons 
voir effet semblable produit pour neutrons. 


246810 i4 6820 24 26 32 36 40 44 48 52 


vallée des masses les nombres protons qui correspondent satur- 
ation des enveloppes des sous-enveloppes protoniques, c’est dire, aux 
nombres magiques que pourrait appeler ‘‘sous-nombres magiques’’. 
Cette figure montre fonction nombre protons largeur 
des éléments pairs. Avec Guggenheimer (6) nous entendons par cette 
grandeur différence entre nombre maximum minimum neutrons 
pouvant liés élément donné. Une telle représentation surtout 
tage bien montrer constance largeur isotopique, d’un 
interval vallée masse qui correspond une enveloppe, une sous- 
enveloppe une somme sous-enveloppes protoniques données. d’autres 
termes, Fig. indique que largeur isotopique varie pas d’une 
enveloppe, c’est dire cours saturation enveloppe, sauf comme 
nous signalé plus haut pour 48, pour cette derniére 
exception étant due neutrons. Naturellement comme pour 
Fig. l’ordre succession des enveloppes adopté dans Fig. est celui 
Maver (7). 


Discontinuités produisant saturation des enveloppes neutroniques 


qui produit pour nombre magique protons voit apparaitre plus 
important décalage tout tableau périodique savoir celui qui lieu pour 
neutrons. Cette nouvelle discontinuité bien visible sur Fig. été re- 
produite détail sur Fig. Cette derniére figure montre qu’a partir 
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unités, limite inférieure unités. d’autres termes, cette dis- 
continuité, non seulement est différente celles qui ont été étudiées 
parce qu’elle produit cette fois pour nombre magique neutrons non 
pour nombre magique protons, mais aussi parce qu’elle s’opére dans 
sens décroissance nombre isotopique 


Les lignes pointillées épaisses Fig. indiquent que serait 
vallée masse réajustement produit pour protons, c’est dire, 
nombre magique protons comptant états. 


84 


schéma Mayer prévoit, aprés nombre magique neutrons, deux 
sous-enveloppes comptant états possibles. Or, trouve que pour 100 
neutrons (82 100) vallée masse semble subir (5) une nou- 
velle discontinuité s’effectuant également vers droite, comme dans cas 
discontinuité précédente, mais moins prononcée. Cette discontinuité ap- 
parait d’ailleurs comme certaine considére les droites pointillées pente 
qui relient les isotopes impair qui brisent avec 100 neu- 
trons. Elle affecterait aussi limite inférieure stabilité 
hafnium avec 100 neutrons, 172, était stable comme cela semble probable 
selon raisonnement Duckworth (4). 


Avec 106 neutrons sous-enveloppe neutronique suivante états est 
saturée; semble encore qu’on puisse établir une correlation entre fait 
déplacement vers droite limite supérieure stabilité qui produit 
pour 106 neutrons qui est visible sur Fig. 
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Enfin, des lignes pointillées fournit que pour 110 neu- 
trons fond vallée des masses déplace vers droite unité. 
nombre coincide avec celui pour lequel saturation sous-enveloppe 
neutronique aprés 106 neutrons, comptant quatre états, est at- 
teinte. est bien connu qu’aucun des noyaux contenant plus que nombre 
magique 126 neutrons n’est stable, aussi semble-t’il priori plus difficile 
suivre les limites extérieures stabilité dans cette région vallée des 
masses les comparer aux limites valables dessous 
naturelle une ligne qui est comparable limite supérieure stabilité des 
noyaux stables, savoir celle qui relie les isotopes Ra226, Th232, 
stables par rapport seulement radioactifs avec une 
trés longue période (1). (Sur Fig. ces noyaux aussi que tous les autres 
noyaux ci-dessous dans région naturelle ont été in- 
diqués sous forme cercles blancs.) Cette ligne trouve une droite 
pente qui prolongée vers bas, recouvre pas avec ligne stabilité 
qui intéresse les noyaux avec 126 mais situe droite cette ligne. 
méme circonstance retrouve pour fond vallée des masses. effet 
Seaborg coll. (8) ont récemment publié une liste des noyaux stables dans 
région naturelle montrant qu’a d’astatine les 
éléments impair ont seulement seul isotope stable. Ces isotopes ont 
été reliés par une ligne pointillée sur Fig. s’apergoit que ces noyaux 
situent nouveau sur une droite pente décalée droite par rapport 
ligne valable dessous plomb, décalage produisant pour 126 neu- 
trons. Finalement liste Seaborg permet considérer, dans région 
naturelle les noyaux stables pauvres neutrons qui sont com- 
parables aux isotopes tombant sur ligne inférieure stabilité des autres 
régions. Ces noyaux, indiqués sur Fig. situent encore sur une droite 
pente subissant elle aussi déplacement d’une unité droite compter 
126. 

III. Discussion 
générale 

L’étude précédente permet faire quelques remarques générales sujet 

des enveloppes nucléaires sur vallée masse. 


fait trés significatif semble étre que saturation des enveloppes, aussi 
bien protoniques que neutroniques, coincide pas dans tous les cas avec 
événement marqué dans vallée; ces exceptions s’observent pour 
région gouvernée par ‘‘régime pente pour 112 (fin des sous- 
enveloppes neutroniques 1h9/2 3p1/2) dans région: pente 
n’est pas correlation avec saturation d’une enveloppe nucléaire. 


Tenant compte ces exceptions semble qu’on puisse avancer 
selon laquelle les régions vallée masse qui correspondent saturation 


a 
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d’enveloppes jouent régions privilégiées d’une dis- 
continuité trouve favorisée supposant toutefois que réajustement 
rapport N/P voisinage ces positions fait sentir. rapport, point 
saturation d’une enveloppe nucléaire correspond aux meilleures conditions 
stabilité requises aucun effet discontinuité cet endroit visible. 
Par contre rapport N/P devient trop inadéquat pour assurer stabilité 
nucléaire cours d’une enveloppe, discontinuité apparait 
avant saturation celle-ci. Cette hypothése est proposée pour 
rendre compte discontinuité produisant pour 46. peut re- 
marquer que cette discontinuité s’est produite cours d’une 
comptant états possibles. saturation était 


produite pour ces deux valeurs. 


qui concerne différence dans des discontinuités qui 

produisent dans premiére partie vallée masse, point saturation 
des enveloppes protoniques, dans deuxiéme partie point saturation 
des enveloppes neutroniques, elle est une étroite corrélation, comme été 
déja remarqué, avec nature régime suivant lequel s’édifient les noyaux. 
effet que ‘‘régime pente favorise l’apparition des discontinuités 
pour les nombres magiques protons tandis que dans 
domaine validité pente toutes les discontinuités (sauf pour 
50) produisent pour les nombres magiques 
neutrons. 


peut aussi faire une remarque rapportant fréquence 
discontinuités produisent avec une plus grande fréquence dans premiére 
partie vallée masse que dans seconde. peut voir preuve que 
‘‘régime pente qui gouverne les noyaux lourds est mode d’édification 
des noyaux qui vérifie davantage les diverses conditions stabilité que ré- 
gime pente fait dans région des noyaux légers moyens. 


qui vient d’étre faite sur correlation entre saturation des en- 
veloppes nucléaires des discontinuités dans vallée masse 
suggére quelques idées nouvelles dans domaine radioactivité naturelle 


Rappelons d’abord que dans région naturelle nous avons 
constater que vallée masse obéit régime ‘‘pente une dis- 
continuité produit pour nombre magique 126 neutrons qui déplace les 
limites extérieures stabilité d’une unité droite par rapport aux limites 
valables pour plusieurs éléments dessous plomb. question qui pose 
tout naturellement est savoir pour combien noyaux avec 126 
régime pente continue sans subir une nouvelle discontinuité qui selon 
toute vraisemblance produit suivant axe horizontal, pour 
nombre donné neutrons (cf. sections IIB III). 
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schéma Mayer (7) prévoit aprés nombre magique 126 nucléons 
une sous-enveloppe comptant états Fig. montre cependant 
qu’a saturation cette sous-enveloppe, pour 138 neu- 
trons, aucun déplacement vallée masse produit. Cette circon- 
stance n’a cependant rien entre dans cadre des exceptions dis- 
cutées dans section saturation d’enveloppe s’accompagne pas 
discontinuité dans vallée masse. D’autre part, examine les 
étendues des intervals durant lesquels les limites extérieures stabilité con- 
tinuent sous forme droite pente sans subir discontinuité arrive 
conclusion que entre les conditions stabilité régime pente 


est d’autant meilleur que les noyaux soumis régime sont plus 


Comme dessous 126 observe une région les limites 
extérieures stabilité poursuivent sans discontinuité pour neutrons 
consécutifs, peut interval continu encore plus étendu 
dessus cet axe, semble improbable avant que 
soit atteint nombre neutrons moins égal 126 146. 


Tenant compte raisonnement qualitatif déduisant schéma 
Mayer que seconde sous-enveloppe aprés nombre magique 
126 neutrons contient états possibles, est conduit envisager nombre 
148 neutrons (126 10) comme pouvant correspondre une nouvelle 
discontinuité dans vallée masse. 


Fig. représente dans plan les noyaux stables autour 
148; ceux prévus par Seaborg coll. (8) ont été suivis points 
rogation. voit que les résultats expérimentaux concernant cette région 
transuranienne vallée masse sont insuffisants pour que puisse 
prononcer d’une maniére certaine sur d’une discontinuité pro- 
duisant pour 148 neutrons. Par contre, est remarquable que les prédictions 


Fic. probables stabilité des noyaux transuraniens. 


144 
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faites pour les noyaux contenant 148 neutrons plus suggérent 
d’une telle discontinuité, qui sur horizontal 148 neutrons, pro- 
duirait d’une unité vers droite serait fait comparable celle qui 
existe pour 126 neutrons. Fig. montre les limites extérieures stabilité 
indiqués suivant les conventions adoptées pour Fig. dans d’une 
telle discontinuité. constate d’abord que seul isotope prédit comme 
stable, 247, place sur nouvelle ligne stabilité maxima non 
sur plus, Seaborg prédit comme plus léger isotope stable Cf, 
246, qui également tombe sur ligne inférieure stabilité dans 
d’une discontinuité produisant pour 148 neutrons. Finalement 
250 qui aurait stable l’ancienne limite supérieure stabilité 
était valable figure pas dans liste precitée. Notons toutefois que Seaborg 
prédit comme stable 244 qui visiblement sur Fig. situe 
dehors nouvelle limite supérieure stabilité. 


Remarquons terminant que nouveau décalage pro- 
duisant pour 148 neutrons suggére revision des prédictions sur les propriétés 
des éléments transcuriens, prédictions faites admettant implicitement 
continuité des limites stabilité dans cette région extréme vallée 
masse. 


particulier, par extrapolation, Seaborg coll. assignent 6.8 6.6 Mev. 
comme énergies des particules des isotopes 246 248 Cf. dans 
thése d’un déplacement vallée des masses pour 148 neutrons les isotopes 
précités placeraient plus prés fond vallée des masses qu’en 
d’un tel déplacement qui est immediatement visible sur Fig. consé- 
quence (ref. 28) les énergies leurs particules seraient plus faibles que 
6.8 6.6 Mev. respectivement. C’est méme probablement détermination 
éxacte des énergies des particules des isotopes transcuriens qui permettra 
premier lieu par comparaison avec les prédictions Seaborg prononcer 
sur validité envisagée ici. 
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NUCLEAR DATA FOR LOW POWER RESEARCH REACTORS* 


The atomic energy authorities Canada, the United Kingdom and the 


United States have agreed the public release certain information low- 
power research reactors, including those nuclear properties uranium im- 
portance the design and operation such reactors. The following are the 
officially released data: 


Thermal Neutron Cross Sections for Uranium 


The following are currently accepted values barns for approximately 
Maxwellian neutron spectrum with most probable neutron velocity 
2200 


Thermal neutron 


cross section for Natural 
Fission 545 3.9 
Capture 100 2.6 3.3 
Scattering 8.2 8.2 8.2 


Natural Uranium Fission Cross Section between 0.7 Mev. and 5.0 Mev. 
The cross section for fission natural uranium neutrons having energies 
between 0.7 Mev. and 5.0 Mev. given the accompanying figure. The 
curve has been drawn fit the most reliable current experimental data. 


Neutrons per Thermal Neutron 


Prompt Neutron Energy Spectrum 
The energy distribution prompt neutrons resulting from the thermal 
neutron fission U-235 given the formula: 

where the neutron energy Mev. the laboratory system. This dis- 
tribution function represents experimental data neutron energies 
Mev. with maximum deviation 15%. 


Fast Fission Effect 

The following are typical values the fast fission contribution the 

reactivity research reactors: 

Ina reactor the type: 2.9% 

Note: and are low-power graphite, uranium reactors 
the Argonne National Laboratory and Harwell, respectively. 
and are low-power heavy water, uranium reactors 
Argonne and Chalk River, respectively. 


N.R.C. No. 2332. 
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Resonance Absorption Integral 
approximate empirical formula for the effective value the resonance 
absorption integral natural uranium 
where the value barns, the integral over the range neutron en- 
ergy from fission energy thermal energy, and where uranium sur- 
face area and uranium mass grams. 

hoped that during later declassification actions appropriate credit may 
given the large number individuals the U.S., U.K., and Canada who 
have contributed knowledge the above nuclear data. 

Values the neutron cross sections for uranium energies other than 
thermal will not released this time except for that information released 
under item above. Similar data other fissionable nuclei will also remain 


FISSION CROSS SECTION OF NATURAL URANIUM FROM 0.7 MEV. TO 5.0 MEV. 


FISSION CROSS SECTION (BARNS) 


0.1 0.2 0.3 0.4 0.5 0.6 0.70.80.910 Ls 20 25 30 35 4.0455.05.5 6.0 


NEUTRON ENERGY (MEV. ) 


Fic. 
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